The Quantum Aspects of Relativistic Fermion Systems 
with Particle Condensation 



OS 

os 

> 
o 

(N 



> 
SO 



Os 



Oh 



X 



1_ 
FT 



S. Ying 

Phystcs Department, Fudan University 
Shanghai 200433, People's Republic of China 

A consistent local approach to the study of interacting relativistic fermion systems with a con- 
densation of bare particles in its ground or vacuum state, which may has a finite matter density, is 
developed. The attention is payed to some of the not so well explored quantum aspects that survive 
the thermodynamic limit. A 4-vector local field, called the primary statistical gauge field, and a 
statistical blocking parameter are introduced for a consistent treatment of the problem. The effects 
of random fluctuations of the fields on local observables are discussed. It is found that quasiparticle 
contributions are not sufficient to saturate local observables. The property of the primary statistical 
gauge field are discussed in some detail. Two models for the strong interaction are then introduced 
and studied using the general framework developed. Four possible phases for these models are found. 
The possibility of spontaneous CP violation and local fermion creation in two of the four phases is 
revealed. The implications of the finding on our understanding of some of the strong interaction 
processes are discussed. 
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I. INTRODUCTION 



Fermions are the fundamental building blocks of the observable universe, which, at certain level, is supposed to 
be understoodable in terms of quantum field theories (QFTs) like QCD, standard model of electroweak interacting, 
etc.. They are, however, less familiar to us theoretically compared to bosons partly due to their lack of classical 
correspondences. An understanding of Nature therefore requires a more direct understanding of the behavior of 
fermions. 

Represented by anticommuting Grassmann numbers in the path integration formulation of a QFT, fermions are 
harder to handle in numerical simulations (like the lattice ones) than bosons. It is therefore desirable to integrate 
the fermionic degrees of freedom out (in the path integration sense) analytically. This turns out to be an easy 
task, at least formally, in most of the cases which deals with a Lagrangian density that is only quadratic in the 
fermion fields. The result is usually a much more complicated effective action for the bosonic fields of the system 
to be functionally integrated by various means like a numerical calculation or simulation, an approximated analytic 
computation, modeling and possibly the mixture of all of them. The fermion loop effects in numerical simulations are 
in some sense less well understood than their bosonic counterparts in a theory at present. The problems become more 
sever in the presence of a finite chemical potential in lattice simulations M. This situation calls for more analytic 
efforts to understand the effects of fermionic quantum fluctuations in an interacting system since it indicates that our 
understanding of the problem is still insufficient. 

The traditional treatment of the finite density problems (at finite temperature) in statistical mechanics is based upon 

the grand canonic ensemble in which the partition function is Z = Tre~P( H ~^ chN ^ with (3 the inversed temperature, 
H the Hamiltonian and N the particle number of the system. A global chemical potential p c h is introduced to select, 
among all possible particle numbers, the corresponding set of particle numbers that are different from each other 
by a finite quantity in the thermodynamic limit. Since only those extensive quantities that are proportional to the 
(infinite) volume of the system are relevant, the above-mentioned differences are irrelevant to bulk thermodynamic 
quantities. This makes the grand canonic ensemble equivalent to the canonic ensemble j^] in which the number of 
particles are kept fixed. The usually calculated quantities, called the apparent particle number here, are expressed as 
N a pp = /3 _1 91n Z/dfi c h is formally identical to what is called the absolute particle number N a b s = J d 3 xTrp(x,t — 
0^ e -P{H-n ch N) i _ y\ can be realized that the identification between N app and N a bs is not mathematically warranted 
in the thermodynamic limit since the particle number N is a macroscopic operator with eigenvalues proportional to 
the volume of the system, we thus expect that the formal equivalence between N app and N a b s and many other physical 
observables so computed to be broken down under certain conditions especially when the quantum fluctuation effects 
are taken into account. In order to characterize the deviation of the quantity N app from N a b s , a dark component for 
each physical local observables in a relativistic QFT is introduced. For example, the dark component of the fermion 
number density operator is defined as Ap — (N a b s ~ N app )/n, with the volume of the system. The questions to be 
assessed are whether or not Ap = and what is the origin of the dark component when Ap ^ 0. 

Such a possibility is important to study because some of the questions posted for the vacuum state of a relativistic 
system governed by certain QFT are different from the ones asked for the condensed matter system in which the 
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quantities under study like the average particle number density is finite with its absolute values playing no direct 
physical role in physical processes and in which the spacetime resolution (energy) of the observation is usually low. 
The vacuum state of a relativistic system, especially the trivial one, is characterized by its nothingness nature, 
namely, all physical observables in the trivial vacuum state are by definition zero. In the non-trivial vacua of the 
system, certain physical quantities develop finite values which need to be evaluated correctly. Those quantities 
like the fermion (baryon) number density and associated energy density should in principle manifest themselves in 
gravitational process at the macroscopic level. In addition, global quantities have no direct physical meaning in a 
classical relativistic system according to principle of relativity. It is expected that the apparent quantities like N app 
are not sufficient ones in the study of the vacuum state of a relativistic system. Rather, one should go back to 
the absolute quantities like N a b s defined above. Therefore, for a better marriage between relativity and quantum 
mechanics, instead of a global chemical potential fj,^, a localized quantity called the primary statistical gauge field 
H a (x) seems to be necessary, which leads to a local theory for the finite density problem. The functional derivative 
of the logarithm of the new partition functional with respect to at certain spacetime location does give the 

absolute particle number density since it is a finite number. 

The principle of locality has far reaching consequences in the development of modern physics. Implied in Maxwell's 
equations for electrodynamics, it motivated the birth of relativity in which it is raised to a principle that governs 
all physical laws in classical physics. Localities in quantum field theories are implemented in most of the theories 
regarded as fundamental like the quantum electrodynamics and quantized non-abelian Yang-Mills gauge theories of 
the standard model. Locality in the quantum field theories, including the fundamental ones, for non- vanishing matter 
and energy density is not, as a matter of fact, fully implemented. Such theories contain inconsistencies, at least at 
the conceptual level, that have to be removed. 

Locality of a symmetry generates a gauge one. For the U{\) invarince corresponding to the conservations of fermion 
number, a new local symmetry called the statistical gauge invariance with the gauge field fJ- a (x) is induced after its 
localization. 

The introduction of a primary statistical gauge field is expected to produce a series of new problems and opportu- 
nities. One of the goals of this paper is to solve these problems and to explore such opportunities so that to develop 
a consistent framework beyond the quasiparticle picture using which the problems related to the vacuum state of a 
relativistic fermion system can be systematically tackled. 

One of the most important problems in understanding a system governed by a QFT that contains interaction is 
to determine the phase structure of its vacuum. The vacuum is a state of the system that has the lowest energy 
that can be different from the trivial one for interacting systems. The non-trivial vacuum of an interacting system 
covered by this study are the ones that contain macroscopic condensation of particles in various form. It has zero 
overlap with the corresponding trivial one when the thermodynamic limit is taken. Such a state can not be reached 
by a perturbative computation, which is local in nature and can only change finite number of particles. A large set 
of non-trivial vacua of interest are expected to be describable in terms of a set of parameters that characterize the 
macroscopic condensation of bare particles. Some of such parameters are called the order parameters since they are 
indicators of a spontaneous breaking down of certain global symmetries of the system. They are accompanied by 
massless Goldstone bosons that generate long range orders which stabilizes the symmetry breaking states. Others, 
together with the order parameters, specify the macroscopic condensation of bare particles in the non-trivial vacuum 
in a more detailed way. For example, the vacuum of the light quark system (up and down quarks) is known to be 
condensed with quark-antiquark pairs. This phenomenon, which induces the spontaneous breaking down of a chiral 
symmetry, is shown to happen in the lattice QCD simulation Q and is supported by experimental facts due to the 
success of the partial conservation of axial vector current (PCAC) relationship and resulting current algebra 

For a system with a large mass gap (compared to the typical excitation mass scale of the system) , the condensation 
of bare particles is unlikely to occur in its vacuum state. So, for the purpose of this paper, I consider light (compared 
to the typical excitation mass scale of the system) fermion systems. They can be approximately represented by a 
massless fermion system. 

An interaction amongst massless fermions and antifermions of the right sign and magnitude is expected to generate 
a non-vanishing number of fermion-antifermion pairs from the bare vacuum. Under certain conditions, the number 
of such pairs can become macroscopic (or proportional to the volume of the system) in the thermodynamic limit. In 
such a case, it is expected to has a phase transition. Such a phase transition was show to happen to the Nambu Jona- 
Lasinio (NJL) model in the quasiparticle approximation (the meaning of which is going to be specified in the following 
sections). It is also shown in Ref. [0J|] that fermion pairs (and antifermion pairs) can condense to lead to a different 
phase (/3 phase) that belongs to the same chiral symmetry breaking chain, namely SU(2)l x SU(2)r — > SU(2)y- 

In general, the condensation of fermion pairs in the vacuum of an interacting fermionic system belongs to one of 
three categories: 1) condensation of correlated fermion and antifermion pairs 2) condensation of correlated fermion 
pairs (and possibly some correlated antifermion pairs 3) condensation of correlated antifermion pairs (and possibly 
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some correlated fermion pairs) . The spin, flavor and color combination of the condensed pairs determines the nature 
of the symmetry breaking channel of the non-trivial vacuum on a finer basis. Since a macroscopic number of fermions 
and antifermions are pumped out of the bare vacuum in the non-trivial vacuum, it is natural to ask what are the 
effects on the its energy density by considering not only the contributions from the long distance or/and time interval 
correlated quasiparticle excitations, but also, at least partially, the contributions from some of the transient and short 
distance quantum fluctuations within the system. To investigate the later effects within a consistent QFT, certain new 
theoretical concepts and tools in describing and interpreting the related physics, to which this work put its emphasis 
on, is proven necessary. 

The paper is divided into three major parts. The first part consists of sections || and VI] , which gives an introduction 



and a summary. In the second part, which consists of sections ||, |l| and IV, a consistent general approach to the 
relativistic fermionic system at both zero and finite density is developed. The third part includes sections ^ and VI , 
in which two 4-fermion interaction models for the strong interaction are introduced and studied using the method 
developed in part two; some novel properties of these models are revealed using the new tools. 

The more detailed arrangement of the paper is given in the following. In section ||, the general framework used to 
handle the fermionic system is discussed. An 8-component "real" representation for the fermion fields is adopted. The 
distinction between the Minkowski and Euclidean spacetime formulation of the problem is emphasized. It is pointed 
out that in the Euclidean spacetime formulation of the problem, additional contributions due to certain quantum 
penetration of field configurations to classically forbidden region, to which quasiparticle approximation in Minkowski 
spacetime can not access, can be included in the effective action. I also motivate the need for a distinction between 



local and global observables in a relativistic QFT. The existence of the dark component is demonstrated. Section [II 
is devoted a tentative local approach to the finite density problem for a Lagrangian density that conserves the fermion 
number. Such a formalism is used to study the question of whether or not a state with non- vanishing fermion number 
density can has a lower energy (density) when the vacuum of the system is in a phase different from the trivial one. 
Spontaneous CP violating stationary points in the Euclidean spacetime is shown to exist in a phase, called the a 
phase of a massless fermion system, with non-vanishing ipip vacuum expectation value. It is argued that such a result 
is at least physically not acceptable. Based on these findings, a consistent theory is developed in section |v| for a 
relativistic fermion system in which a phase transition with particle condensation has occurred. A statistical blocking 
parameter e is introduced. The question of the statistical gauge invariance of the system due to the original global 
U(l) symmetry related to fermion number conservation is addressed. To demonstrate the relevance of the theory, two 
models for strong interaction are introduced in section Their phase structures are then studied. The vacuum of 
both of them has three different phases. One is the trivial (bare) vacuum, which is called the O phase; the second 
is the above mentioned a phase with fermion and anti-fermion pairs condensation; the third kind of phases, called 
the f3 phase and ui phase respectively, are phases that spontaneously break the global U(l) symmetry related to the 
fermion number because of a condensation of fermion pairs and antifermion pairs. These phases are further analyzed 
in section |Vl| by using the formalism developed in section IV . A spontaneous CP violation is found to be present in 
the (3 and u> phases. Also, the spontaneous creation of matter, which is relevant to Cosmology, is s hown to appear 



naturally in the (3 and uj phases. Finally, the main results are summarized and discussed in section VII which also 
contains an outlook. 

II. GENERAL DISCUSSIONS 
A. Fermion representation 

The fermion field is represented by an 8-component "real" spinor 'F . It is essential for a consistent formulation of 
a relativistic finite density field theory in a functional (or path integration) approach developed in this work Q (see 
also Appendix |X|) and for other aspects of the QFTs || . In this representation, 'F satisfies the reality condition 



where the fig matrix is 



tf(a;) = y T (-x)n , (2.1) 



<>„ ()>C- ( ^ C Q 1 ] (2.2) 



with C the charge conjugation operator. The matrices O-y is one of the three Pauli matrices Oi,2,3 acting on the upper 
and lower 4 components of $ . 
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For the case of the number of flavors n / for the fermions is less than three, the symmetry transformation of the 
8 component \P can be made to be t he s ame as the 4 component representation ip by imposing a slightly different 



"reality' condition on <!> given by Eq. 2.1 with SIq replaced by 



Here p— liinf — 1 and p = 1T2 if n/ = 2. More detailed properties of this representation for the case of rif = 2 are 
discussed in Refs. They will not b e repeated here. 



If n f > 3, the re ality condition Eq. 2J5 is no longer valid. In such a case, the general "reality" condition is given 



by the original Eq. 2.1 and the representation of VP under flavor SU(rif) transformation generated by 



T _ f t a 3 If t a is symmetric , , 

a 1 t a If t a is antisymmetric ^ ' ' 

have to be adopted. Here a = 1, 2, . . . , n?j — 1 and i a is the generator for the symmetry transformation in the 
4 component representation of ip. The set of matrices T a and t a (a — 1,2, ...,n/) belong to equivalent adjoint 
representation of flavor SU(nf) transformation. 

It should be emphasized that the intention of introducing an 8-component representation for the fermions is different 
from the one related to the doubling of degrees of freedom in the closed time path approach (h],[11|] to non-equilibrium 
(also equilibrium) problems or thermal field dynamical |l2| approach to equilibrium problems in two aspects: l) this 
work is devoted to the study of zero temperature physics where, in the language of the closed time path approach, all 
the dynamical fields considered in this work lies on a single time axis running from negative infinity to positive infinity 
rather than on two time axis that form a closed loop 2) there is no doubling of degrees of freedom for the fermions 



here since the constraint Eq. 2.1 is systematically implemented in the formalism. In case the present formalism is 
to be extended to finite temperature H or to the non-equilibrium situations, a doubling of the components of the 
8-component "real" representation has to be made in addition. 

B. Minkowski spacetime formulation 

The generating functional can be written as 



e w[J,v, v ] 



J [[Dl/ilfllfle'/^^^ 1 ^^ 8 '^^^* 4 ^), (2.5) 



where J = { J\, J2, . . . , J n } are a collection of external fields coupled to the corresponding boson fields / = 
/2, ■ ■ ■ j /«.}) £b[/] * s ^ ne Lagrangian density for the boson fields /, and ij, fj are external Grassmann fields 
coupled to the fermion fields ^. Here /j can be either real or complex. W[J, rj, rj\ generates the Green functions 
of the fermion or the boson fields. The possible gauge fixing conditions, which can be implemented by multiplying 
1~L D[fi] a set of S functions or by introducing ghost fields |13| , are unimportant to this work and are suppressed in 
the sequel. 

The bosonic part of the Lagrangian density will not be specified in this investigation. This allows the results of this 
work to be useful in a wide class of problems. For example, in case of QCD, the boson fields are 8 gluon fields B^(x) 
fa = 1, 2, . . . , 8) with the full Lagrangian density in the 8 component representation for quarks provided in Appendix 

In the study of the vacuum properties, fermion degrees of freedom can be eliminated first 



J ]^D[f i }e^ LnDet ' <olS ^ 1U]+ ^ SFU], ' +i I d4x ( CB ^+T, k Jfc ^-). (2.6) 



Here "Dei" denotes functional determinant. The proper vertex generating functional for the boson fields is then 

r[/] = w[j,o,o]-iJ2 J *f k - ( 2J ) 

k 

The stationary configuration / determined by the equation 

ffil=0 (i=l,2,...,n) (2.8) 

OJi 
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with fi (i = 1,2, ... ,n) spacetime independent determines the phase structure of the vacuum. In general, T[f] is 
difficult to compute directly. It is useful to define an effective action S e f / for the boson fields as 

Seffif] = -i^LnDetjoiS^tf] + i^LnDet l0 iSp l [Q] + J d 4 x£ B [f] 

-i^SpLnS^if] + i^SpLnSp^Q] + J d 4 xC B [f], (2.9) 



whe re S p denotes the functional trace and a constant (infinite) term independent of the boson fields are subtracted. 



Eq. |2.6| becomes 

e W[J,rj, v ] = J JJ D [f.yS e ffU]+hvS P U]v+i J d4 *Y, k J *f«. (2.10) 

i 

Starting from S e f / [/] , either systematic improvement beyond the Gaussian approximation or numerical simulations 
like lattice computation can be made. A more detailed discussion of the local quantum fluctuations around the mean 
field is given in Appendix |b|. A formal relation between r[/] and S e ff[f) is developed in Ref. |Q and discussed 
in Appendix When the fluctuation in / is only treated at an one loop level, the vertex functional T[f] is (see 
Appendix ^) 

= iSefAf) - \spLnDG- l [f], (2.11) 

where D is the bare propagator for / and 5 2 S e ff/SfSf is symbolically denoted as G _1 [/] . Under such an approxi- 
mation, the solution to the equation 

%^ + >G[/]^l=0 (; = l,2,...,n) (2.12) 
of j 2 5 f J 

determines the vacuum phase structure of the system. It will not be further discussed in this paper since the results 



of this paper depend only on some of the bosonic fields {fi} as a solution to Eq. 2.8 being different from zero. 

The effective action S e ff[f] can be expressed in terms of the spectra of the operator joiSp which is Hcrmitian. 
The eigenequation of interest is 

7 ^ f 1 [/]*a = A[/]vI/ a (2.13) 
with ty\ the eigenvector. In the time independent case, S e ff[f] in terms of A is 

^Af] = -%^^§ ]+ J^ Blf] , ,2,4, 

where T is the temporal dimension of the system (T — > oo), p° represents the energy of the eigenvector ^\ and £ is 
a collection of other quantum numbers that completely determines a single eigenvector The order in which p° 
integration and £ summation is carried out is important in general since they both are divergent before the subtraction. 
The symbol 

is understood as that neither the sum over £ nor the p° integration is not done first but they are done in a covariant 
way. The order between them depends on situations which are discussed in the following and in Appendix |A|. 

Due to the logarithmic function, the integrand in the complex p° plane is multivalued, the physical contour C with 
Feynman-Mathews-Salam causal structure [0 for the p° integration is shown in Fig. [l]. The conventional computation 
of the energy density of the vacuum with a non-covariant cutoff, which is called the quasiparticle approximation, can 
be represented by a distortion of contour to curve I of Fig. ^| It corresponds to a summation of the energies of 
individual stationary quasiparticle orbits in the negative energy Dirac sea with a fixed (time-independent) background 
configuration for /. 
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C. Euclidean spacetime formulation 



The path integration representation of the generating functional W[J,Tj, rj\ on the right hand side of Eq. 2.5 
contains ambiguities associated with the non-specification of the initial and final field configurations. For the transition 
amplitude between a given pair of initial and finial field configurations, the contributing intermediate states can in 
principle be different from the vacuum state interested here. The lowest energy configuration corresponding to the 
vacuum is automatically projected out in an Euclidean spacetime computation with a sufficiently large Euclidean 
time for a large set of proper initial and final fie ld c onfigurations. 

In addition, the Minkowski effective action Eq. |2.9| has extrema (or is stable) only at configurations of the boson fields 
that are of classical nature. An important class of configurations corresponding to pure quantum mechanical effects, 
namely the tunneling effects through potential barrier or penetration into the classically forbidden field configurations 
are expected to be missing in a steepest-descent or Gaussian approximation. The Euclidean action obtained by 
replacing time variable t by —it (with i 2 = —1) can be stable at either the time independent configurations, which are 
the same ones as in the Minkowski approach, or the configurations that correspond to the tunneling or penetration 
that are quantum mechanical in origin. It is important to realize that the later stable configurations are always 
absent in the set of stable configurations of the Minkowski effective action despite they constitute an important 
contribution to the energy density of the system. The Euclidean spacetime formulation has been used to find important 
stable field configurations called instantons in non-Abelian gauge theories. They correspond to stable finite action 
field configurations due to the tunneling between gauge field configurations of different winding numbers. Other 
applications of finding the tunneling effects by using the Euclidean spacetime formulation and their connection to the 
WKB method in quantum mechanics can be found in e.g. Ref. jTsj. It will not be elaborated here. The effects that 
are interested in this study are the ones that survive the thermodynamic limit and are thus non-finite action effects 
despite they can be decomposed into a collection of random finite action ones (see Appendix [b]). 

The generating functional for the Green functions in the Euclidean spacetime formulation can be expressed as 

e W[J,v, v ] = Jl[ D[f l ]e S "ff [f]+ ^ s " [f]r,+ I di ' x T. k J »f* (2.15) 

i 

which also serves as a definition of the Euclidean effective action S^j [/] . The Euclidean propagator S§ [/] is found 
by using the rules discussed in the following. 

For a fermion system, the simple replacement t — > —it is not sufficient to obtain a consistent Euclidean spacetime 
formalism from the corresponding Minkowski one due to the fact that a fermion is not a scalar particle. More 
sophisticated set of transformations are needed. Formally, an Euclidean effective action for fermions can be obtained 
by making a continuous change of the metric Eq ]. The result of change for a Dirac particle can be summarized by 

«& = {-.--.-}. 7tf = *7 5 , li = ~Ho, if =H (2.16) 



and the effective action Eq. 2.9 in the Euclidean spacetime formulation becomes 



S eff[f) = \ {SpLn [ l75 (5|[/])- 1 ] - SpLn [*76(5f [0])" 1 ] } + J <fix E L B [f E \, (2.17) 



where the Euclidean propagator Sp[f] is also obtained by using substitution rules listed in Eqs. 2.16 . In terms of the 
eigenvalue of the hermitian operator i^iSplf]) -1 that satisfies 

*75^|[/])- 1 *a = A[/]* a . (2.18) 



The Euclidean correspondence of Eq. 2.14 is then 



where the superscript "£"' on top of a quantity in Euclidean spacetime is suppressed in there and in the following 
whenever no confusion is thought to occur. It can be demonstrated that Eq. 2.19| can be obtained from Eq. [2.14 by 
a distortion of the p° integration contour to curve II shown in Fig. ||. 
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D. Local and global observables 



The physical observables in a classical theory consistent with relativity are local ones like the charge density, energy 
density etc.. Unlike in the non-relativistic world, the global observables like the total charge and the energy of the 
system are not directly measurable. In a given frame, however, the global observables can be defined as a spatial 
integration of the local observables on an equal time hypersurface in the Minkowski spacetime. The meaning of the 
global observables can only be defined operationally. The particular value of a global observable has to be determined 
by the following procedure. First a synchronization of the clocks of a group of observers on each spacetime point 
should be carried out, then let each observer measure the corresponding density at his/her location in spacetime at 
the same time and finally sum (integrate) each observers finding to obtain the value of the global observable. 

The global observables in quantum theories, including the relativistic QFTs, are regarded as physical observables. 
Let us consider certain density observable p(x) in a QFT. The corresponding global observable Q is defined simply as 

Q= I d 3 xp(x) (2.20) 
Js 

with S the equal-time hypersurface in certain reference frame. Then the matrix elements of p(x) and Q, namely, 
(/ |jo(x)| i) and (/ \Q\ i) define the corresponding observables. 

Since both the local observable p(x) and the classically not defined global one Q are defined in a relativistic 
QFT, one can naturally ask the following question, namely, is there any difference between O — (Q,\Q\Q.) and 
O' = J d 3 x{£l \p(x) \ ft) measured in a state fi? It can be shown that this is a relevant question for a relativistic QFT. 

A measurement of the total charge or charge density of a system normally consists of supplying an external global 
field or local field coupled to the corresponding observables in a known way and then measure the response of the 
system like the force that the external field exerts on the system or the amount of increase of some proper defined 
"potential" of the system, which allows the observer to deduce the total charge or charge density measured. Let us 
consider the measurements, in the above sense, of the global and the local observables in the vacuum state of the 
system expressed symbolically as 

Q = lim(0. 7 \Q\ 0j) = \im(0j | f d 3 xp(x)\ 0j), (2.21) 

O' = I d 3 x Jim (0 5j(x) \p(x)\Q sj{x) ), (2.22) 

where J is a global external field coupled to Q and 5j(x) is a local external field taking non- vanishing value only at 
x that couples to p(x) and | 0j) and | 0sj(x)) are the corresponding vacuum states. The first one Oq corresponds to a 

measurement of Q directly; the second one corresponds to the integration of a (infinite) set of measurements of p(x) 
on a space-like hypersurface. These two measurements can in principle be different in a system governed by a QFT 
because of the random local quantum fluctuations of the fields that are not suppressed in the thermodynamic limit . 

Such a potential difference guarantees the possible existence of the dark component in a QFT. The fact that 
localized random quantum fluctuations are not suppressed in an interacting QFT, especially in some of the non-trivial 
phases of the system, are discussed in more details in Appendix |b] where it is also shown that the conventionally 
used quasiparticle picture in many-body theory and QFT is not sufficient to saturate the local observables due to the 
existence of the dark component. Such a deviation from the conventional physical picture based upon quasiparticles 
gets less and less significant as the resolution of our observation respect to spacetime gets lower and lower compared 
to the typical size of the localized random fluctuations f a of the system. In such a case these fluctuations are more 
and more suppressed and the contribution of the dark component becomes smaller and smaller resulting in I) an 
emergency of a quasiparticle dominated picture for the system and 2) the validity of the results obtained based upon 
a global chemical potential p c h in the grand canonic ensemble in low resolution (energy) observations. 

III. A NAIVE LOCAL QUANTUM FINITE DENSITY FIELD THEORY AND ITS PROBLEMS 

A. A tentative formalism for the local finite density fermionic field theory 

In order to develop a theoretical framework consistent with the requirement of locality, a new quantity called 
primary statistical gauge field p a {x) is introduced in the following. The motivation for its introduction is discussed 
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in the introduction and in the above sections. The other reason that it is treated as a local variable is due to the fact 
that total fermion number has limited meaning in a relativistic system that are generated in the past not infinite long 
ago. The total number of fermions in such a system is not an observable since there are regions outside the horizon 
that are classically non-detectable even in principle due to the constraint of causality. 



1. The asymptotic grand canonic ensemble 

The generating functional corresponding to Eq. |2.5| for a fermionic system with finite density can be formally 
written as 

e W[J,rj,w] = J TJ^[J i ] jD [^] e «/^(^- 1 [/]*+Mor+£ B [/] + -..) ) (3 J-) 

i 

where ". . ." denotes the source terms and the fermion number current is 

j»(x) = \^{x)^0^{x). (3.2) 

In the Minkowski spacetime, quantity on the left hand side of the above equation is the transition amplitude of the 
system from properly weighted initial field configurations 4>i at t = — oo to final field configurations (f>f at t = +00, 
namely, 

e w[J,n,n,n] = ^ W[0 /> 0iK0 /j i = +oo|0 4> t = -oo) J (3.3) 

where W[(f>f, 4>i\ is the weight functional discussed in the following. 

These initial and final fields are considered free fields. The interaction terms are adiabatically switched on at certain 
large negative time — T and switched off at certain large positive time T . Such a technical manipulation does not 
affect the actual local physical observables like the energy density, fermion number density, etc. at time t that is far 
from both — T and T due to locality. 

Since the particle content for free fields at t — ±00 is clear, which allows a straight forward statistical interpretation 
in terms of number of particles in each single particle state of the system. It is natural to assume that the initial 
(final) state are in the grand canonic ensemble with the weight W[0/, <fii] for the summation determined by the factor 
lim^oo cxp[— f3(E — /iiV)][J where E is the total energy and N is the total fermion number of the (free) system, fx 
agrees with the time component of the spacetime independent part of the primary statistical gauge field fx a (x). Such 
an ensemble is called the asymptotic grand canonic ensemble here. 

The asymptotic grand canonic ensemble differs from the grand canonic ensemble for it allows for a local approach to 
the relativistic finite density problem and the existence of the dark component discussed in the above sections and in 
Appendix [b|. Its predictions, however, approaches that of the grand canonic ensemble in the non-relativistic situations 
in which the spacetime resolution or energy in a measurement is low. This point is demonstrated in Appendix |b] at 
the mean field level. 

Two points must be pointed out. The first one is that although the time interval 2T between the (adiabatic) 
switching on and off of the interaction terms are let to go to infinity in the final result, the thermodynamic limit, in 
which the spacetime box that contains the system approaches infinity, is taken first. Therefore T is still infinitesimal 
compared to the temporal size of the system. The second one is that, as a result, for the conserved quantities, like the 
total fermion number, the extreme value of it picked out Q in the asymptotic grand canonic ensemble with a fixed 
value of (i is unmodified. The later point is elaborated in the following. 



1 We are interested in the zero temperature case here. The form of the weight functional for free fields at finite temperature 
can be found using the method given in, e.g., Refs. jHjjjnJ, which involves two time axis: one runs from —00 to +00 on real t 
axis; the other lies below it on t he complex t plan. In the zero temperature limit, only the contributions from the real t axis is 
nonzero, which leads to Eq. 3.1. 



9 



2. The fermion number density in the asymptotic grand canonic ensemble 



Before continuing the development of the formalism, it is important to reveal a property of the asymptotic grand 
canonic ensemble related to the U(l) symmetry corresponding to the fermion number conservation. The conservation 
of j a (x) due to the £/(!) symmetry can be derived from the Noether's theorem at the classical level, namely, 



d^(x) = 0. 



(3.4) 



Those interaction Lagrangian densities for which Eq. 3.4 remains true at the quantum level are considered despite the 
fact that this symmetry is not explicitly related to a gauge symmetry for which a superselection sector in the Hilbert 
space of the system exists. The reason to consider such a class of models is because for a quark system, the fermion 
number is identical to the baryon number, which is conserved to a high precision due to the lack of any convincing 
evidence of the proton decay i n ob servation at the present. 
For an uniform system, Eq. 



3.4 implies that 



dgi 







(i = l,2,...) 



(3.5) 



in the asymptotic grand canonic ensemble with {gi} representing a set of interaction coupling constants and the 
derivative taken by keeping /i Q unchanged. For an uniform system, the mean local fermion density is a function both 
of the coupling constants {gi} and /i Q , which is now spacetime independent. Eq. 3.5 implies that the mean local 
fermion density of the system ~p is only a function of fi a , it is independent of the interaction coupling constants {gi}. 
This property allows us to find out the relationship between p and [i a easily by considering a non-interacting system. 
It is discussed in Appendix [A]. The result for a massless fermion with n / flavors and n c colors is 



P : 



UfUc 3 
37T2 M ' 



(3.6) 



where /i = \J [)? ■ The simplicity of Eq. |3.6| in the asymptotic grand canonic ensemble is due to the fact that the 
interaction Lagrangian density conserves fermion number. It implies that for an interacting massless fermion system, 
local quantity p is non-zero as long as fi is non-zero even when a phase transition in its vacuum state that generates 
a finite gap for the lowest excitation of the system has happened. Such a qualitative behavior is required following 
the discussion given in Appendix 0. Eq. 3.6 is however an exact relation in the asymptotic grand canonic ensemble. 



After an exploration of the implications of the U(l) symmetry of the class of models under consideration in the 
asymptotic grand canonic ensemble, we are in a position to further devel op t he formalism for the inves tiga tion of 
the density fluctuations of the vacuum after a phase transition. Since Eq. 3.1 can be obtained from Eq. 2.5 by the 
replacement 



ifi + /i0 3 



the effective action S e f / of the system is changed to 



K°Af;tA 




V f — l 



d 4 xC B [f] 



with the eigenvalues A obtained from the following equation 

7°^ 1 [/;m]*a = A[/;m]*> 



(3.7) 



(3.8) 



(3.9) 



where S'^ 1 [/;/Li] is derived from the corresponding one in Eq. 2.13| by making the substitution Eq. 3.7. Since the 
question interested in this study is related to a comparison of the energies of states with different fermion densities when 
the interaction is present, I consider the following effective action obtained from S e ft by a Legendre transformation 



S e ff - S, 



eff 



,4 -c 
d Xfl a J 



(3.10) 



It is a canonic functional of j a with /i a implicitly depending on it. For an uniform system like the vacuum, the 
second logarithmic term in Eq. 3.8 is calculated in such a way as not to over-co unt the already known (and included) 
quantum fluctuations of the free fields. The result is given in Appendix [A] (Eq. A7), it is 
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^J c ^j^m=J dx ^- £) (3 - n) 



in the rest frame of the density. Together with Eqs. |3.6| , Eq. |3. 10 for an uniform system takes the form 
~ dp A p o e [/;/i] n f n c f 4 4 f 4 

^ = - Z 2S^ ln V^lM-l^i +Jd X C B [f}. (3.12) 



The effective potential to be minimized for an uniform system is then defined by 

^e// = - lim Seff/VzT 
V3T— >oo 

1 ^ dp \,o^[f;iA n f n c 4 



where volume V3 = L 3 with L the spatial dimension of the system. 

B. Euclidean Instability of the a Phase and Its CP Problem 

The a phase for a massless fermion system can be shown to realize by using the Nambu Jona-Lasinio (NJL) model 
p7| | . In the 8-dimensional "real" representation for the fermion spinor, the NJL model with 2 flavors (rif — 2) and 3 
colors (n c — 3) are given in appendix |j. The NJL model is studied extensively in the literature by assuming /x" = 
in the vacuu m so far. It's possible that such a plausible assumption is in fact incorrect. 

With Eq. |3. 13 , the question of whether or not the a phase is stable against fluctuations in \i a can be studied. In 



the a phase, Eq. p. 13] takes the following form 

(2tt)*-V- p\) \ p 2 J + 4G a + 2tt 



where p+ = (p° + fi°,p + fi) and pt = (p° — fi°,p — fi). 

1. Quasiparticle approximation and its problem 

In the quasiparticle approximation, the p° integration contour is the one shown in Fig. |^. The localized quantum 
fluct uations of the order parameter a are not included following such a contour. With the help of Appendix [A], Eq. 



3.14| can be shown to have the following form 

3 4 1 4 6 



V e ff = - ~U + 



2tt 2 ^ 4G 



^2 



£ 3 d\p\\p\ 2 (VIpP + ^-W), (3.15) 



where fc^ = — a)^ /i 2 — u 2 with 9{x) the step function and A3 is the cutoff in the 3-momentum that defines 
the model. It is a monotonic increasing function of fi. So the stationary point for V e ff is at the position /i = 0, as 
expected. 

There is however an inconsistency related to the fermion number in the quasiparticle approximation to the a phase. 



Eq. 3.6 holds for any phase of a model Lagrangian density for which the fermion number is conserved. On the other 



hand, the average fermion number for an uniform system in the quasiparticle approximation is 

«4%- ff )(/' 2 - ff2 ) 3/2 ( 3 - 16 ) 



which differs from Eq. 3^ and the results of Appendix [b] even qualitatively. The reason, as is mentioned generally 



in section |II], is because the quasiparticle approximation is related to a Gaussian approximation to the generating 
functional of the model in the Minkowski spacetime formulation. Such an approximation can only take into account of 
the contributions from propagating modes or quasiparticles. The quasiparticles are not fundamental building blocks, 
namely, the bare particles, of the system; they are composite objects representing coherent excitations of infinite pairs 
of bare particles and anti-particles. In the phase that quasiparticles propagate, the excitations correspond to the bare 
particles of the system damp in spacetime when produced. They are not, however, absent in a time interval that is 
sufficiently short. They form an important component, which is called the dark component (see Appendix |b[), that 
contributes to the fermion number density and other relevant local physical quantities. The contribution of such a dark 
component is expected to be taken into account, at least partially, if we formulate our semi-classical approximation 
in the Euclidean spacetime to sample the important contributions of pure quantum mechanical configurations. 
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2. Euclidean stationary points and the CP problem 



In the Euclidean spacetime, Eq. 3.14 can be evaluated by using the p° integration cont our sh own in Fig. S together 



with a covariant cutoff. This operation corresponds to the replacement rule given by Eq. 2.16 plus /i° — > , which 
is equivalent to the change p° — » ip° and fi° — > fi° in the expression for the Minkowski effective action after the tracing 
over spin, isospin and color degrees of freedom is taken. The time componen t of ii a is kept unchanged since it is 



regarded as an external field and is related to the fermion density given by Eq. 3J3. The result is 

Veff/A = ^Jo ^ Jo dxVl - x ln tf-JPF+APyW + T6^ (T (3 - 17) 

(7 = cr/A, jl = fi/A, ao = GoA 2 /47r and with A the covariant cutoff in the Euclidean momentum space. 

In the a phase with non-zero cr, the dependence of V4///A 4 on /i is plotted in Fig. [| The minima of V e ff is not 
located at fi = for any finite a but some finite < < cr. This is physically not acceptable. 

There are two problems related to such a vacuum with non- vanishing fermion density if the dominate phase of the 
physical strong interaction vacuum is considered to be in the a phase. 

The first one is related to the strong CP problem. Since in such a background field as /i", the two CP conjugate 
eigenstates of a neutral particle has different energies; it implies the occupation number for one eigenstate can be 
larger than the other in physical processes at sufficiently low energy. This in turn will cause much too large CP 
violation phenomena not observed in nature. Let us consider the only system, namely the K° / K system, in which a 
CP violation was observed. A mechanism for the explanation of the CP violation in the K {) / K system was proposed 

long ago |pj| . Its basic idea is to assume a \i a like potential through out the space that differentiates K° and K . 
However, such a fjP is estimated to have a value of order 10 -8 eV. It is much smaller than the average value found 
here, which is of the order of 10 2 MeV if A is taken to be 1 GeV. 

The second one is related to the dark matter problem. If there exists an uniform fi° ~ 10 2 MeV field in the universe 
at the present, the dark component of the baryon number density would be of order of the nuclear matter density, 
which is much larger than what it is expected. Had such a scenario been correct, the universe would not has lived 
until today. 

Facing these two serious problems, a solution is needed to be looked for. There are at least two alternatives. The 
first one, which is likely to be correct, is that there are something missing in the computation procedure used so far. 
The second one is that our notion that the dominate phase of the universe is in the a phase is wrong. The later 
alternative is unlikely to be correct since there is a large body of empirical facts that support such a notion. 

With this consideration in mind, I turn next to an investigation of the first alternative. 

C. A Fock space inspection of the vacuum structure of the a phase and the blocking effects 

It is known that in the a phase of the NJL model treated in the mean field approximation, the vacuum can be 
related to the bare one by an unitary transformation before the thermodynamic limit L 3 — > oo is taken. It can be 
explicitly written as 

| vac) = JJe^P'^lO) (3.18) 
p . h 

with p the three momentum, h the helicity label and operator O expressed in term of the creation operators (a ph , 
b^-ph) an d annihilation operators (a p /t, &- p a) of the bare fermions as 



0(p,h) = -e p al h bl ph -b- ph a ph , (3.19) 



where cos# p = |p|/_E p and E p = \/\p\ 2 + cr 2 . The annihilation operators for the quasiparticles a p h and the anti- 
quasiparticle j3 p h of the system are related to the bare ones through a Bogoliubov transformation 

a ph = cos(i6> p )a p/! . - sin(i0 p )&L ph) (3.20) 
Pph = sin(^e p )al ph + cos(ifl p )6p h . (3.21) 
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From Eqs. 3.18, 3.19 two properties can be seen: 1) the a phase vacuum is a superposition of states in the Fock 
space of bare fermions with increasing number of fermion and anti-fermion pairs and 2) after the thermodynamic 
limit L? — > oo is taken the overlap between the bare vacuum | 0) and the a phase vacuum | vac) tends to zero. 
This implies 1) the occupation of fermions and anti- fermions in the a phase vacuum can have blocking effects on 
the creation operation of bare particle states and 2) the true a phase vacuum can not be reached by perturbative 
iterations starting from some preassumed state of the system without introducing some kind of macroscopic variables 
into the path integration formalism. 

IV. A CONSISTENT LOCAL QUANTUM FINITE DENSITY FIELD THEORY 
A. The statistical blocking parameter 

1. The motivating Fock space study 

I develop a theory that takes into account both the fermion-antifermion pair condensation and the resulting block- 
ing effects in this section. As mentioned in the previous section, the bare vacuum with zero pairs of fermion and 
antifcrmions is not a suitable initial state that starts the path integration computation since it has zero overlap with 
the true vacuum state of the system after the phase transition. Instead, the starting state is of the following kind 

with n the number of the fermion-antifermion pairs and £ other quantum numbers that completely specify the state. 
The diagonal matrix clement of the evolution operator of the system can be expressed in terms of path integration 
over the dynamical fields of the system 

(0 o ,t = +oo|^ o ,f = -oo)= lim (0o|e~^ ( '' _tO |to) 

t f oo 
ti — > — oo 

= E E C i C (* «< * = +°o | n£, *=-«>) + ..., (4.2) 

n «' 

where H is the total Hamiltonian of the system and . . . represents the off diagonal contributions to the transition 
amplitude between states with different pairs of fermion and antifermion. For a discussion of the eigenstate of the 
total Hamiltonian like the vacuum state, the off diagonal contributions with macroscopically different initial and final 
fermion-antifermion pairs are expected to be effectively absent in the final result after the thermodynamic limit. 

The transition amplitude (n£', tf | n£, U) with the external fields present is then written in terms of path integration, 
namely, 

«',< = +oc|<,i = -oo) =N j ^[^TJ^f/.jei/^Kl^F'I/l^+^I/l+^+W+E.^A), (4.3) 

where N is a constant. The formal manipulations, which express the above functional integration over fermion degrees 
of freedom by an effective action S e ft, remain mostly unchanged, except S e ff depends now on a new statistical 
parameter e 

t = +oo | n£, t = -oo) = N' J Y[ D { fi yS cff lM+^s F lf]v^fd* x J2 k J »h . (4.4) 



For a stationary situation, the effective action S e ff[f, /x, e] is given by Eq. 2.14. The constraint that both of the initial 
and the final states considered are configurations with both of the bare fermion and antifermion states below energy 
e filled is implemented by a distortion of the contour for p° integra tion from the one in Fig. to that of Fig. § 



In the thermodynamic limit of L 3 — > 00, the sum over n in Eq. 4.2 can be replaced by an integration over e, the 
statistical blocking parameter, namely 



E 



deM(e) (4.5) 
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with M(e) the integration measure of the transformation. Eq. 4.2 becomes 

(4> , t = +00 \(f> ,t= -00) = / deM(e)^2CICp{n4',t = +00 \n£,t = -00} 

= J de J Y[D[fi\e iSc " [f ^' e]+ ^ SF[f]v+l I dix (E k J "fk-Vo^,e)) 



(4.6) 



with exp[— i J d 4 xVo(^, e)] the leading piece of the measure for the e integration in the thermodynamic limit. When the 
volume L 3 of the system becomes increasingly large, the integrand of the e integration becomes increasingly sharp at 
the extrema positions of the argument of the exponential in the above equation. This is due to the fact that e couples 
to macroscopic variables that are proportional to the volume; it has no quantum fluctuation in the thermodynamic 
limit. Due to this reason, the detailed form of the measure M(e) but its leading piece in Eq. 4.5 is irrelevant in the 
thermodynamic limit. So the weighted sum 



W[<fo,M<fo,* = +oo| 



6 , t = -00) = 



(4.7) 



with e taking the value of one of the extrema of the argument of the exponential. The conventional method, which 
may turns out to be not sufficient in symmetry breaking phases, corresponds to a special case, namely, the e = one. 
The relevance of introducing e will be discussed in the following sections. It can be seen that a possible finite e in the 
final result is perturbatively non-reachable. 



2. The determination ofVo(n,e) 



Eq. 4.7 tells us that the generating functional of an interacting fermionic system should be written as 

e W[J,rj.,ri.^.e] = J JJ D ^ e iS eff lf,^e] + irjS F [f}r,+i J d i x{Y Jk Jkfk-V {p.t))_ 

i 

It is normalized by the condition 

W[0,0,0,Q,0] = 0. 

For a stationary situation, S e ff corresponding to Eq. 3.10| can be expressed (see also Eq. 2.14) as 



+ / d 4 x[C B [f}-V (^e)}, 



(4.8) 



(4.9) 



(4.10) 



where C represents the p° integration contour shown in Fig. |^ and Co re prese nts the p° integration contour given by 
Fig. |l} Vo(/U, e) satisfies V^(/u, 0) = J d i x[i a j a so that it agrees with Eq. 3.10 in this special case. From Appendix [A|, 
it is shown that the left hand side of the following equation is finite for an uniform system, namely, 

J d 4 x [(m+P(+) + M-P(-) - Mp) - (e(+) + e ( _) - e) 



The basic assumption of the theory is then the following choice for Vo(/^,e), namely, 

V (fi, e) = / d A x + M-P(-) - MPj 



(4.11) 



(4.12) 
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which completely specifies Eq. 4.7. For an uniform system, the corresponding effective potential in a Hartree-Fock 

(4.13) 



approximation to be minimized is 
V eff = lim 



dp X p a^[f;fj] n f n c , 

o_ ln \ . rn ..1 + /i_2 U 



vwio"2F 3 ^c S 2tt 111 A p o )S [0,/x] 4vr 2 v/ ' 



2e 4 + 12 M V) 



with C denoting the p° integration contour chosen. 

In order to preserve the causal structure of the original Minkowski p° integration contour, the Euclidean effective 
action is obtained by distorting the Minkowski contour given in Fig. || to the one labeled "II" in Fig. |. 



B. More on the primary statistical gauge field 



1. Statistical gauge invariance, physical states and conservation of fermion number 



In the process of introducing the primary statistical gauge field /i Q , the original global U(l) symmetry corresponding 
to the fermion number conservation is replaced, in a certain sense, by a local symmetry. This local symmetry originates 
from the fact that the eigenvalues A p o^[/, fj], which satisfies the eigenequation 



7 ^V,M]*A = A[/,/z]*A, 
is invariant under the following gauge transformation 

*A(a;) -» e l ^° 3 ^ x {x) 
fi a (x) -+ p a {x)-d a 4>{x) 



(4.14) 



(4.15) 
(4.16) 



with 4>{x) an arbitrary function of the spacetime that decreases to zero sufficiently fast at the spacetime infinity. 
The introduction of a local field ^ a {x) is expected to introduce infinite extra degrees of freedom, which should be 
eliminated in certain way. Albeit the full effective action give n b y Eq. 4.10, which depends on p a p a > 0, is not 
invariant under the gauge transformations given by Eqs. |Hj and pLl~6| , the primary statistical gauge invariance of 
the quantum fluctuation or the connected part of S e ff requires further investigation. 

Let us find the connection of the primary statistical gauge invariance to the conservation of fermion number by 
quantize the system governed by the Lagrangian density 



£ = C + fi a Jc 



(4.17) 



in Eq. 4.8. Here C is the original Lagrangian density before introducing p a and j a is the fermion number current 
density given by Eq. 3.2. The fermion field ^ and the boson fields {/,} are quantized as usual ||; they shall not be 
repeated here. 



What is needed to be found here is the conjugate variable 7r" of [i a . For that purpose, Eq. 4.17 can be treated as 
the Hamiltonian density, namely 



= £. 

Before the quantization, it follows from the Hamiltonian dynamics that 

dn (p) 



<9o7T ui 



dpi 
dir m 



: Ji 



0, 



where i,j = 1, 2, 3 labels the spatial components of 4-vectors. 

The quantization of /ij is then implemented by the Dirac quantization condition 

[TT ui (x,t),fij(xf,t}] = -i5^ (x - x')Sij . 



(4.18) 



(4.19) 

(4.20) 
(4.21) 



(4.22) 



The statistical "electric field" it U i commutes with all elementary fields in the Lagrangian density but the one listed 
above at equal-time. Note that all quantities with a hat "A" on top denote operators in the following. 
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After the quantization, the set of gauge transformation, in which <fi(x) is independent of time, is represented by 



* -> f/[0]Wt[0] = e^° 3 * (4.23) 

/I Q -> [f[0]/2 Q [/t[0] = m« - 9 Q (4.24) 

with 

= £ r ird M'+ v ^> (4.25) 

where the space integration is at any specific time at which a transformation of an operator is considered. The 
operator algebra between dynamical fields is then represented in a Hilbert space. Since the quantity p a is introduced 
into the original theory, it is expected that this Hilbert space contains states that are not physical or are redundant. 
The physical states are selected within the full Hilbert space by requiring that they satisfy 

(Phys' \U[(f>]\ Phys) = Ze~^ (4.26) 

under the time independent gauge transformation discussed above, with the common phase factor $ restrict ed to 



those functions that are independe nt of time (it is explained in the following). Taking the time derivative of Eq. 4.26 
and using dynamical equation Eq. |4.19| , one obtains 



(Phys' | (d / 



p + V • j ) | Phys) = (Phys' Phys) = (4.27) 



which is the conservation of fermion number current in physical processes. 

The time independent and spatially localized gauge transformation considered is non-trivial one. It selects amongst 
those states in the extended Hilbert space the physical ones. This can be understood if one consider the commutation 
relation between p + V • n u and the Hamiltonian of the system. 

p + V = ij t @+ V • 7?„) = dj" (4.28) 

with H the total Hamiltonian of the system. It is zero due to the conservation of fermion number current. It means 
that p + V • tt u is independent of time when the matrix elements between physical states are taken. The states in 
the extended Hilbert space can be divided into subspaces labeled by a complex (time independent) function of the 
spatial coordinates according to the matrix elements of p + V • n u between themselves. For those eigenstates of the 
Hamiltonian of the system, it can be written asQ 

(<p\ \(p + V • 7T„)| <4) = 6 EiE] N l3 , (4.29) 

with <j the space dependent complex function, \(p5) a state in the physical space that has energy Ef., See 1 taking 
zero or unity value if E =/= E' or E' = E (assu ming that E i s discrete before thermodynamic limit is taken) and 



Nij independent of spacetime coordinates^. Eqs. 4.28 and 4.2£ mean that the physical states are the ones that have 



vanishing matrix elements on the commutator of the Hamiltonian and jo + V • n u . Therefore they are expressible by a 
superselection sector in the extended Hilbert space defined and labeled by a complex function <j of spatial coordinates. 
Such a definition of the physical states for the primary statistical gauge theory is less restrictive than the one used 
in dynamical gauge theory like QED [^lj where due to the existence of the dynamical part for the gauge fields at the 
tree level, the physical states are restricted to the subspace with <r = only. In fact, for the (3 and w phases discussed 
in the following, in which the U(l) symmetry corresponding to the fermion number conservation is spontaneously 
broken down, c can not be zero due to the fact that before taking into account of the dynamical gauge fields (that of 
the photon) , the massless Goldstone boson corresponding to the spontaneous symmetry breaking has to be considered 
as a physical excitation. But if q is chosen to be zero, the massless Goldstone boson belongs to unphysical states [p2[ . 
Such a situation actually opens up the possibility for the spontaneous CP violation to be discussed i n the following. 



The choices made for the physical states is a constraint invariant under time evolution due to Eq. 4.28. It shows 
that definition for physical states remains true at all times and no transition to unphysical states and between the 
superselection sectors is possible in physical processes. 



2 It is the quantum version of the "classical" constraint equation p + V • 7r u = <;. 

3 Eq. 4.28 may not be restrictive enough. Somewhat more restrictive constraint can be suggested. It consists of decomposing 
7r u into superposition of holomorphic (p + V ■ TV n )r-) and antiholomorphic (jo + V ■ 7f u )( + )Components 



p + V ■ n u into superposition of holomorphic (p + V ■ n u 
states are those ones that are eigenstates of (p + V ■ tt u 



201. The physical 



(-) 



with a common eigenvalue c. 
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2. The statistical gauge degrees of freedom and the question of long range order 



The primary statistical gauge field fj, a is non-dynamical at the tree level. At the quantum level a dynamics for \i a 
is generated due to the fermio n qua ntum fluctuations. The relevant effective action for the primary statistical gauge 
field can be obtained from Eq. 4. 10| by a "Legendre transformation" to a form without the contribution of Vq, namely 



S ef f[f,n,e] 




dp Q A p o^[/, /i] 



«2tt A p o !? [0, M ] 



£ £ / 



(4.30) 



which is a canonic functional of /^". The quadratic term for slow varying /i Q (in spacetime or at long distances) 
generated from the fermion determinant is of the form 



S eff = | / ^xd A X^' a {x)^{x, X ')^{x') + 



1 



n f n c 



a 2 a 

a xe fi , 



(4.31) 



where fi' a = n a —~p, a with ~p a shifted \x a and the last term is from the corresponding one in Eq. 4.30. The first term is 
generated from the fermion determinant. If the electromagnetic interaction between the fermions are not considered, 
TT al3 (x,x') is given by 



nf(x,x') = i(0\Tj a (x)f(x')\0) 

= Z^(d 2 x g^ - 9°dg)5(x - x 1 ) + i(0 \j a (x)\ 0)(0 \f(x')\ 0} 



(4.32) 



in the normal phase; and, 

nf(x,x') = i(0\Tr(x)f(x')\0) 



fj 



gag0 

^X ^x 
Q2 



S(x-x') + i(0 |r(x)|0)(0|/(x')|0) 



(4.33) 



in the phase where the U(l) symmetry corresponding to the fermion number conservation is spontaneously broken 
down since there exists a massless pole in the matrix element of j a (see Ref. for a more detailed discussion). Here 



Z^> and 11^ are functions of x and x' 



In the normal phase, the effective action for slow varying \x' a is 



°eff ~ 



d 4 x 



U«F V + \ (»<0 |r| 0)(0 |/| 0) + 2g«^e 2 ) p! a p!, 



(4.34) 



with = dV' 3 -9V° 

since the eigenvalues A p o j are invariant under the gauge transformation given by Eqs. |4.15| 
and 4.16. In the phase where U(l) is spontaneously broken down and before considering electromagnetic interaction, 



S[f f = ~ Id'x 



{ [i(0 \j a \ 0)<0 |/| 0) + g ati (n<"> + 2^e 2 ) 



(4.35) 



In both of the situations with slow varying fx' a , Z^> and II^) are approximately co nstan ts. 

If the electromagnetic interaction is considered, then there is a f a pf a ^ term in Eq. 4.35 and the II^) term is absent, 
even in a spontaneous U(l) symmetry breaking phase. This is because when the electromagnetic interaction between 

the fermions are considered, 7r Q ^ can be decomposed into a connected and disconnected part or ir$ +i(0 jj^l 0)(0 \ \ 0) 
with the connected part given by 



(4.36) 



and 



( c ) i - \ v 

= (l 9v,v - quq»)n 



(c) 



(4.37) 
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where Gq V is the propagator of the bare photon, — > in the g M — > limit even when Kq} v contains a massless 
pole (in q 2 ). Whatever the case, it can be easily seen that the full two point proper vertex for p! a is non- vanishing in 
the — > limit if either e / or (0 |^| 0) ^ or both in the vacuum state of the system. Therefore a long range 
order for fi' is possible only if both e = and fi = in the vacuum. The spatial component fj,' of ii' a is short ranged 
in the phase where e ^ and /i = 0. It is however long ranged in the phase where e = 0. 



3. Topological configurations 



The discussion given above shows that after the introduction of a primary statistical gauge field fi a , the repre- 
sentation Hilbert space of the operator algebra is extended. Such an extended Hilbert spac e inc ludes not only the 
physical states, but also non-physical ones. The physical states are the ones that satisfy Eq. 4.29 , which is expected 
to be satisfied in all consistent computations of physical observable where fermion number current conservation is 
preserved. 

The extension of the representing Hilbert space gives us additional leverage to project out the collective excitations 
and configurations not easily discernible in the conventional approach. 



Due to the statistical gauge invariance under the local transformation given by Eqs. |4.15| and |4.16| , the physically 
non-trivial configurations of the system depend, in the path integration sense, only on the flux density or statistical 
"magnetic field" b defined as 



b = V x fi. 



(4.38) 



Therefore the general form of the generating functional Eq. 4.8 can be written more precisely by imposing certain 
gauge fixing condition or as 



,w 7 



(4.39) 



where J(b) is the integration measure. Eq. 4.39 is equivalent to Eq. 4.8. It is however useful for us to find out 
collective stationary configurations that are not easily found in a common approach to the effective action. 
For a configuration in which b is non-vanishing only in a localized region, quantization of flux appear, namely 



dS ■ b = d> dl ■ = 2nir, 
s J 9S 



(n = 0,±l,±2,...), 



(4.40) 



where E is the surface that contains the b and line integration is around the edge <9E of S. The quantization results, 
as it is well known, by imposing the uniqueness condition on the eigenfunctions IP \ [/, /j] . 



4- A new macroscopic parameter and long range order 

Using the primary statistical gauge field /i Q , a new macroscopic parameter that characterizes the vacuum state of 
the system can be introduced. It is defined as 

dv = e i $ dl ' fJ> , (4.41) 

where S is a large 2-dimensional surface area in space at certain time and the line integration <f is along the edge of 
the area E. 

It is known that the vacuum expectation value of Os provides another one of the macroscopic parameters for a 
more detailed characterization of the phase structure of the system. For example, in a disordered system in which the 
correlation between the complex phase of the fermions at different space points becomes short ranged, condensation of 



vortices or monopoles of the type of configurations with non-vanishing n in Eq. 4.4C can derive a Kosterlitz-Thouless 
pij type of phase transition. 
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V. TWO MODELS FOR STRONG INTERACTION AND THEIR VACUUM PHASE STRUCTURE 



The fundamental theory for the strong interaction is considered to be QCD with its Lagrangian density given in 
Appendix |c]. The current available method of studying QCD from first principle is the lattice QCD simulation. Albeit 
great progresses are made, the lattice computation are limited by the small lattice sizes and by the limitation in the 
computer power. Model approaches, which are simpler than the full QCD calculation and has large overlap with it in 
the low and intermediate energy regions can be and have been used to obtain much of the physical pictures that are 
supposed to happen in the system governed by QCD Lagrangian density. Since the mass of the light quarks have values 
much smaller than the typical scale of the hadronic spectrum of order 1 GeV, certain subset of the behaviors of the 
light quark system can be simulated by an interacting massless fermion systems that possesses the basic symmetries 
of the QCD Lagrangian density. It is expected that we can learn some of the important possible behaviors of the 
light quark system by using models due to our (relatively) increased theoretical analytic power. The spontaneous 
chiral SU (2) j, x SU (2)# symmetry breaking down in hadronic systems was historically discussed before the birth of 
the conception of quark and QCD Lagrangian density. This phenomenon is only latter justified by the QCD (lattice) 
calculation. 

I discuss in this paper the physical properties of strong interaction vacuum related to the fluctuation of baryon 
number using two model Lagrangian densities that possess the chiral SU(2)l x SU{2)r symmetry of massless QCD 
and has 3 colors (n c = 3). For a full investigation of the possible phases of the vacuum of a relativistic massless 
fermion system, these models are so chosen that they allow not only the quark-antiquark condensation that is widely 
discussed in the literature but also the rarely studied phases that are induced by a condensation of quark-quark (or 
antiquark-antiquark) pairs. One of these possibilities is studied in detail in Ref. |7],^j . 

In order to simplify our discussion, I consider two model Lagrangian densities that are half bosonized. Both of them 
can be viewed as been the descendents of two four quark interaction models after introducing the auxiliary fields in 
a Fierz invariant way @]. 

This section, which mainly serves the purpose of introducing the models, contains the determination of their vacuum 
phase structure using the conventional approaches. A more detailed study using the refined method developed in 
sections III and IV is given in the next section. 



A. Model I and the uj phase 

The first model is defined by the following Lagrangian density 

A = ~* [*? - <r - «? ■ r 7 5 3 - 7 5 Ax c O (+) - 7 5 ^ c X c O ( _)] * - ^-(a 2 + tt 2 ) + ^XcX' \ (5.1) 

where a, ir, x c and \ c are auxiliary fields with (x c )^ = — x c an( i ^o, G3/ are coupling constants of the model. A c and 
A c (c— 1,2, 3) act on the color space of the quark; they are 

A C C1C2 = -e CClC2 A C , C1C , = e CClC2 (5.2) 

with e abc (a,b — 1,2,3) the total antisymmetric Levi-Civita tensor. Here the quark spinor is represented by the 
8-dimensional Dirac spinor and 0(±) are raising and lowering operators respectively in the upper and lowering 4 
components of "J. 

The effective action is given by Eq. 2.14 with the auxiliary fields independent of spacetime, the effective potential 
has the following form 

1 '»» ^rE ln T^)+ii CT2 -^r^ ( 5 - 3 ) 



2l,t-, oL 3 T^ \(°) 4G 2G 3 , 



where A ra and An correspond to the eigenvalues of the two Hermitian operators defined in Eq. 2.13 with and without 
the auxiliary fields shifted respectively. Since the auxiliary fields do not depend on spacetime, the eigenvalues X n and 
A n ° can be labeled by the 4-momentum of the corresponding eigenstates ^x n . The result is 



eff 



* f d4 P foV-i X *(P) , X '^P) \ , 1 2 1 
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where \i(p), X'iip) are eigenvalues of states with color different, the same as \ c respectively and factors 8, 4 correspond 
to their degeneracy. It has the following explicit form 



V e 



eff 



d 4 p 



In 



4Gn 



2G,, ' 



(5.5) 



where x 2 = -X C X C - 

In the Minkowski spacetime formulation, the contour for the p° integration is shown in Fig. [I]. One of the most 
commonly used ones, which is called quasiparticle path, is shown in Fig. |2[ It can be shown that the resulting 
effective potential is the change of the energy density of the quasiparticles, which is obtained from summing over the 
quasiparticle energies determined by the poles of Sf [/] , relative to that of the bare particle in the truncated Dirac 
sea. It does not has a covariant form due to the fact that a non-covariant cutoff in the 3-momentum p space has to be 
introduced. The Euclidean path shown in Fig. || can be used to obtain a covariant expression for the effective action. 
The Euclidean covariant cutoff scheme does not result in deriving physically undesirable results in other covariant 
approaches |p5| . As discussed in section 2.2, it can also be used to include quantum effects (see Appendix |b|) that are 
beyond the quasiparticle approximation. The differences between these two paths are however not very important for 
the purpose of this subsection. Since different kind of cutoffs are used in these two approaches, a comparison between 
them is difficult. Nevertheless, I shall adopt the Euclidean path shown in Fig. g. An Euclidean path is however 
important for discussions to be followed. The resulting expression is 



d 4 p 

WT 4 



In 



9 i 2 \ 2 



p 



p 



p 



4G 



2G,< ' 



(5.6) 



where A is the Euclidean cutoff introduced to define the model. A numerical evaluation shows that the minima of 
V e ff(<7 7 x) is located on either the a axis (x = 0) or the x axis ( CT = 0). Explicit expression for V e ff(a, 0) and V e ff (0, %) 
are found to be 



Ve// (a,0)=3/(— ) + 



1 G 1 



A 2 ' 167ra A 2 



(5.7) 
(5.8) 



where the dimensionless effective potential v e ff is defined by V e // = A 4 v e ff, ao — GoA 2 /47r and ay = GyA 2 /8tt 
with 



m = A 

87T Z 



-x + In [ 1 + - ] x 2 - ln(l + x) 



(5.9) 



The values of a 2 and x 2 a t the minima of Eqs. 5.7 and 5.8 determine the vacuum of the system in the one loop 



Hartree-Fock approximation for the fermions (i.e. Eq. |2.1l| without the second bosonic one loop term). The phase 
structure of the model is shown in the ag-ay plane in Fig. kU. Three kinds of phases for the vacuum are possible. The 
first phase, which is called the O phase, is the bare vacuum. The second phase, or the a phase, has non-vanishing 
vacuum expectation value of \E'\E'; the chiral SU(2)l X SU(2)n symmetry is spontaneously broken down to a SU(2)v 
flavor symmetry. The third phase, called the lo phase of the vacuum, has non-vanishing diquark and antidiquark 
condensation characterized by a non- vanishing x 2 ', chiral symmetry is unbroken in this phase. 

The phase transitions across the boundary between the O and the a phases (ao = 7r/12 and ay < tt/8) and the 
one between the O and the u> phases (ao < 7r/12 and ay = tt/8) are of second order. The phase transition between 
the a and the u> phases (ao > 7r/12 and ay > tt/8) is of first order. The Meissner effects for the electromagnetic field 
are expected in the u phase. The basic physics of it is discussed in more detail in Refs. |f7|,p|Jl8| for model II of the 
following. I shall relegate such a discussion for this model to other work. 

The Minkowski propagator for the quarks in the lo phase can be found by an inversion of the operator in 
the action, namely. 



S F = i[ip- i 5 A cX c O {+) - 7 f U c x c O(-)] 
In the momentum space, it can be expressed explicitly as 



(5.10) 
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Sf(p) 



{i> 7 5 Ax c O (+) - 7 5 A c x c O { - } , ) -jrf-* 



7 



For quark type I 
For quark type II 



(5.11) 



where a quark of "type I" is the one that has color different from \ c ( or X c ) an d a quark of " type II" is the one 
that has the same color as x c ( or conjugate to that of \ c )- It can be seen that a quark of type II has no gap for its 



excitation and an excitation of a quark of type I has a finite gap \J \ 



B. Model II and the (3 phase 



The second model Lagrangian density is 



J0.2 = -V P 



if -a -in- T7 5 3 + (^YVA + <5£ • tV A) (+) 



1 



{AC 



f-f 5 A c + <V • f 7 M c ) O 



'(-) 



„ -(a z + n z ) - ^U, 
4G 2G 3 



(5.12) 



where cr, ff, <£°, C „, J° and (5 CM are auxiliary fields. It is symmetric under the chiral SU(2)l x SU(2)r group 
transformation. This model is discussed in detail in Refs. [0J|]. It is written down here for references in the later 
sections. The Euclidean effective potential for this model used for a determination of the vacuum phase structure in 
a Hartree-Fock approximation is found to be 



V 



eff _ 




) 2 + -(l + -)(l 



-4(1 + 4>G# 



(5.13) 



A numerical evaluation of Eq. 5.13 shows that the absolute minimum of V e ^(a 2 , cf) 2 ) is located at either a 2 ^ and 
2 = or a 2 = and 2 ^ in the spontaneous symmetry breaking phases. The phase diagram for this model is 
presented in Fig. [5]. Three kinds of phases for the vacuum are possible. The first phase, which is identical to the O 
phase discussed above, is the bare vacuum. The second phase, which is the same as the a phase introduced above, 
has non- vanishing vacuum expectation value of Vl/Vl/; the chiral SU(2)l x SU(2)r symmetry is spontaneously broken 
down to a SU(2)y flavor symmetry in this phase. The third phase is labeled as the /3 phase. It has non-vanishing 
diquark and antidiquark condensation characterized by a non- vanishing (f> 2 ; the chiral symmetry is spontaneously 
broken down to a flavor symmetry, in the same way as in the a phase. 

The phase transition is of second order across O phase and the a phase boundary (a Q — n/12 and a 3 < 7r/4). It is 
first order phase transition across the a and (3 phase boundary (pta > tt/12 and a 3 > 7r/4). There is a second order 
phase transition across the O and /3 phase boundary (ao < n/12 and a 3 = 7r/4). 

In the f3 phase, the propagators for the quarks are found M to be 



Sf(p) = { 



where 4> 2 = <^ c <^° = -<^0 Mt 
non- vanishing auxiliary fields 



■l 

p 



For quark of type I 
For quark of type II 



(5.14) 



In order to simplify the computation, a special choice for the complex phase of the 
'■ and <f> is made, namely, 4> cu is chosen to be </> = ~<t>u- 
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VI. THE VACUUM STRUCTURE OF THE a , /3 AND w PHASES 



The possible phases of two model Lagrangian densities are studied in the previous section using the conventional 
approach to effective potential in the Hartree-Fock approximation. A more detailed characterization of these vacua 
and their properties is given in the following using the general framework developed in section IV. 



A. The a phase 



The effective potential for t he a phase in the full theory given in section [iy| can be obtained by computing the right 



hand side (r.h.s.) of Eq. 3.14 using the p° integration contour C in Fig. || namely 



d A p , _f 1 a 2 \ a 2 \ 1 2 , 3 



V ei M e ) = 6 V c (# 1 H 1 "^J [ 1 -pZ) + 4Go a+ 5?" ■ (6 - 1} 

The effective potential has only a trivial minimum located at n = e = when the integration contour is chosen 
to be the quasiparticle one shown in Fig. ^. This contour, which selects only the quasiparticle contributions, violate 



the conservation of baryon number explicitly; it is demonstrated in section III. It is necessary to turn the contour C 
into the complex plane to find the Euclidean stationary points of the field configurations in a way that preserves the 
causal structure of Fig. |(| Such a contour is also displayed in Fig. 0. 

A numerical study shows that the minima of V e f / for a given value of non-zero a is located either at /i ^ and 
e = or at /i — and e ^ 0. Fig. || shows V B ff along three different directions in the // — e plane. The absolute 
minimum of V e ff is located at e = e$ ^ and /x = 0. This is the result I regard as natural since it avoids the not 
observed CP violation associated with the a phase found in section It also agrees with the physical picture for the 
a phase, which is condensed with correlating quark-antiquark pairs. 

With a non-vanishing e, Let us first assess the nature of the primary statistical gauge field excitation. The effective 
action for the primary statistical gauge field fi' a at long distances is given by (see Eq. 4.34) 



4 Vtt 



(6.2) 



Since e 2 ^ in the a phase, the /i' a excitation is massive (or short ranged) in the static limit and is now stable against 
quantum fluctuations. This agrees also with our observations since no corresponding long range force and large CP 
violation are observed at the present-time condition. 

Due to the presence of a finite e, the response of the system to external (or internal) excitations is different from the 
familiar one we have learned from the conventional approaches. When Dirac's view for the bare vacuum of fermions is 
taken, namely the bare vacuum corresponds to a state in which the negative energy states are filled and the positive 
states empty, the vacuum where e ^ is the state in which all single particle states with energy E < — e and < E < e 
are filled whereas other single particle states empty. Since for an uniform system, the value of the mass a for the 
quasiparticle is always larger than e in the models considered, the presence of a finite e in the vacuum of the a phase 
appears to has no effects if the quasiparticle can propagate long enough without suffering from further scattering. 
The presence of e only provides a virtual possibility for an uniform system that the local fluctuations of the fields can 
feel. So, it has genuine physical effects on local observables even in uniform systems according to the discussion in 
Appendix |B]. 



For a non-uniform system in which the energy of the operator in Eq. 3.9 could be smaller than e, the presence of e 
may has a real effect on the dynamical processes. For example, in a chiral soliton in which the energy of the lowest 
orbits for the valence fermions moves with the size of the soliton, the presence of e will limit the range of change of 
the soliton's possible size which gives an extra stability of such solitons. 

To put the above argumentation in a more concrete context, let us consider a situation in which the lowest energy 
valence fermion state lies within the region — e < E < for one size and shape, then a change in its size or shape 
that moves E upward can be continued freely only until E = since the < E < e states are filled and the next 
available state is the one with energy E — e, which can only be reached by a discontinuous change in the size and 
shape of the soliton. If the nucleon can be regarded as a chiral soliton, this mechanism can prevent it from dissolving 
inside a nucleus if the lowest energy valence (constituent) quark states inside the nucleon lies between (— e,0). Other 
implications of a nonvanishing e are worth to be studied in future works. 

Perhaps other interesting implications of a non- vanishing e are on the particle production and dissipation processes 
in non-equilibrium situations like the heavy ion collision. 
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B. The j3 and lj phases 



According to the physical picture discussed so far, the baryon number content of the j3 and to phases ought to be 
different from the a phase due to the fact that, instead of quark and antiquark pairs, quark pairs and antiquark pairs 
(or diquark and antidiquark) are condensed in the vacuum. This is reflected in the fact that in these two vacua, 
the expectation value of \ c (together with its conjugate field x c ) and (together with its conjugate field 4> cfl ) are 
non- vanishing. 

In the (3 and u> phases, where a = 0, the effective potentials for model I and II are found, using Eq. 4.13, to be of 
the following forms 



V, 
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For Model I 



For Model II 



(6.3) 



(6.4) 



Numerical evaluations show that the local minima of the above effective potentials are located either at fj, ^ and 
e = or a = and £ ^ 0. The corresponding V e f / in the fi-e plane along three different directions are plotted in Figs. 
H and [To] respectively. The absolute minima of both effective potential are located at fi ^ and e = 0. That is, in 
the [3 and lu phases where quark pair or diquark con dense, the vacua of the systems are the ones with finite density 
of baryons with the baryon density given by Eq. 3.6. 



Such a phase also spontaneously violates the CP invariance of the system's original Lagrangian density due to the 
existence of a non-vanishing CP odd order parameter fi a in these phases. In addition, a pattern in which baryonic 
matter and antimatter are separated in space for the (3 and u> phases of the vacuum are energetically favored. The 
superselection sector in the Hilbert space in which the physical states stays in for s uch a CP violating phase can then 
be determined. Assuming that the statistical "electric" field tv u is finite, then Eq. i ,2S| tells us that c = p when the 
spacetime approaches infinity. Due to the translational invariance of the vacuum state, it is natural to require that 
<; = ~p at all spacetime position. In this way the physical states labeled by (, in the CP violating phase of the system 
are determined. 

One of the interesting properties of the (3 and u> phases is that these phases are expected to have off diagonal 
long range order ODLRO p6| due to the spontaneous breaking down of the U(l) symmetry corresponding to baryon 
number conservation. Macroscopic quantum phenomena manifest in those phases possessing ODLROQ. The quantum 
nature of these phases leads to behaviors of the system not expected from our daily experiences. Some of these 
behaviors are observed in the superfluid state of 4 He. Could it allows a quantum mechanical jump (transition) from 
an initially zero baryon density (3 or oj phases of the vacuum state to the lowest energy state of the system inside 
relatively large regions: some of them contain baryonic matter and others contain antibaryonic matter? Such a kind 
of transition is forbidden in classical picture since it violates the the relativistic causality. It is allowed in the quantum 
measurement processes according to the standard interpretation of quantum mechanics, in which a collapse of the 
wave function of the system occurs after a measurement. It remains to be understood in the future more detailed 
researches. If it is permitted, then each of these regions can has a finite size at the moment of the transition or, in 
another word, each of them can has a size large than its event horizon. 

This property of the (3 and uj phases could provide a mechanism for the baryogenesis in the early universe in a 
matter-antimatter symmetric universe p8[ . It was shown to be impossible for such a baryogenesis mechanism to 
be compatible with observations if the process of baryon-antibaryon separation is classical []29|| . The discovery of 
the possible (3 and to phases for the strong interaction vacuum in this study might provide a theoretical basis for a 
reconsideration of the idea of matter antimatter symmetric universe. Some of the other more detailed consequences 
of this picture, which is beyond scope of this paper, is worthy of exploring. 



4 ODLRO is absent in the normal phases of matter where classical picture is supposed to emerge for macroscopic systems 
due to decoherence characterized by a diagonalization of the effective density matrix of the system interested during the time 
evolution (see, for example, Refs. |27|). 
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The //q degrees of freedom is non-propagating in the (3 and m ph ases due to a non-vanishing (uniform) baryon 
number density is present in these two phases. According to Eq. 4.35| 



r/ 4 .; 



(6.5) 



It shows that the /j? fluctuation is damped by a non-oscillating Gaussian factor with the width proportional to the 
inverse square of the baryon number density in these phases (i.e 
phases is however long ranged. This can also be realized through an inspection of Eq 



(?). The spatial component of ji a in the (3 and oj 



4.3£ 



The question of whether or not there is a condensation of statistical mono-poles in these two phases, which is 
discussed in a general term in subsection IV B| should be studied in future more detailed works. 



C. The excitations of the primary statistical gauge field 



In the above discussion it can be seen that in the a , (3 and oj ph ases, either e or fi (it is equivalent to (0 \j \ 0}) is 
nonvanishing. From the discussion presented in the subsection IV B , it can be concluded that there is no long range 
order for the time component of the primary statistical gauge field [i a at long distances or /i° corresponds to at most 
a massive excitation in the non-trivial vacuum discussed in this paper. 

The situation for the spatial component of fi a is different between the a and (3 or uj phases. In the a phase, the 
excitation related to [i a is short ranged. In the (3 and w phases, the excitations related to fj, are long ranged. These 
excitations can therefore generate a statistical "magnetic force" between different particles within the (3 or oj phases 
of the vacuum. The consequences of such a statistical "magnetic force" on the evolution of the system is worthy of 
studying. 



VII. DISCUSSION AND OUTLOOK 



It is found that at least three interwinding new theoretical elements are necessary to be brought into a consistent 
treatment of the problem. The first one is the general existence of the so called dark component in an interacting 
system originated from the transient and short distance quantum fluctuations of the system, which is measured by 
the difference between the absolute value and the apparent value of some conserved quantities like the baryon number 
density, energy density, etc. of the system. The second one is related to the recognition of the existence of the so 
called blocking effects in the non-trivial phases of a system. The third one is related to the necessity of introducing a 
primary statistical gauge field coupled to the fermion (baryon) number current density of the system. By introducing 
these three elements into the formulation of the problem, the door to go beyond the physical pictures limited by the 
approximated concept of quasiparticles is open, which allows us to explore new physical possibilities. 

A systematic path integration formalism for the investigation of the quantum aspects of an interacting fermion 
system (or sector) sampled by Euclidean spacetime stationary configurations in which a condensation of fermion pairs 
(fermion pairs, antifermion pairs, and fermion-antifermion pairs) is present is developed based upon the asymptotic 
grand canonic ensemble. Two statistical parameters, namely the primary statistical gauge field [i a and the statistical 
blocking parameter e are introduced to allow a finer characterization of the vacuum structure of the system. In 
addition, it is shown that the asymptotic grand canonic ensemble reduces to the grand canonic ensemble as the 
spacetime resolution of observation is sufficiently lowered. Such a behavior is a necessary condition for the usefulness 
of describing the macroscopic properties of an interacting system in terms of particles in certain domain of energy 
and for the smooth approach to the well established results in non-relativistic condensed matter systems at low 
energies. Combined with the Euclidean approach to the effective action, some of the quantum effects that survive 
the thermodynamic limit can be included. The present approach, which uses quantum field theoretical language, is 
consistent with thermodynamics and can be extended to finite temperature case Q . 

Firstly, the dark component for local observables like the fermion number density does exist in interacting theories 
in which the direct association of the field theoretic definition of fermion number density with the number of "free 
particles" per unit volume becomes obscure especially when a phase transition inside the system has occurred. This 
conclusion is also applicable to other local observables like the energy density, which may have implications on the 
dark matter problem in Cosmology since it implies that the apparent matterless space at the macroscopic level is 
capable of revealing itself of matter effects in low energy gravitational processes when [i is below the baryonic particle 
production threshold even after the energy density for the [i = state is subtracted. This is because gravitational fields 
couple locally to the source matter fields which contain the random quantum fluctuation generated dark component. 
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Further researches in that direction in the context of understanding the cosmological baryogenesis and dark matter 
problem is an interesting direction to be explored. 

In addition, the picture that a nucleon is made of three valence quarks (quasiparticles) need to be modified when 
there is additional close-by virtual phases for the hadronic vacuum state that has slightly higher energy density than 
the actual one. The implications of such a finding can be explored in observables related to a nucleon. Some of them 
are studied in Refs. pjl8|, other related problems concerning a nucleon, like the understanding of the origin of the 
Gottfried sum rule violation in deep inelastic scattering confirmed in a recent measurement |nj , the small-x behavior 
of nucleon structure functions in deep inelastic scattering [^2],p3[ , etc.. 

Secondly, it is important to take into account the fermionic blocking effects due to the presence of a macroscopic 
population of bare particles in the nontrivial vacuum phases of a system. The blocking effects are generally included 
in the theory by introducing the statistical blacking parameter e, which is non-zero for a system's certain vacuum 
phases. 

The effects of the blocking can not be progressively generated using perturbative expansion starting from a field 
configuration with e = 0. In the a phase of the strong interaction vacuum discussed here, e = configurations are 
inconsistent ones since it is known that the a phase of the strong interaction vacuum is macroscopically populated 
with the current quarks and antiquarks, which changes the available states for an current quark due to Pauli principle. 
The discovery of the blocking effects has hitherto unnoticed implications related to, e.g., the stability of a nucleon in a 
nucleus and nuclear matter, the mechanism for particle production in a heavy ion collision, new ways of (quasi)particle 
dissipation in a strong interacting system, etc.. 

Thirdly, the statistical gauge degrees of freedom of the system represents certain collective mode of the system that 
has a dynamics of its own. There are two components for the primary statistical gauge field: the first one is the 
classical configurations which serves as a background field; the second one is the local fluctuations of it around the 
classical configurations. 

The classical configurations, which is a spacetime independent background [i a in the case studied here, play the 
role of the chemical potential in the the conventional non-relativistic approach. It determines the asymptotic grand 
canonic ensemble of the system. It can also has non-trivial topological configurations corresponding to different 
quantized statistical "magnetic flux" which can determine the phase structure of the system on a finer basis. Once 
present, the statistical "magnetic field" affects the dynamical evolution of the system that can result in 1) material 
pattern formation and 2) providing the seed for the galactic magnetic field in the early universe. Whether or not such 
an idea is actually relevant to comprehend what happened in the the early universe can be studied in further works. 

The local fluctuations of it around the background extended configuration represent the corresponding dynamical 
excitations of the system. The necessary condition for the existence of long range statistical gauge correlation in 
various possible phases of the system is discussed. 

The statistical gauge degrees of freedom are also introduced in condensed matter physics in the context of the 
half filled Hubbard model, which serves as one of the prefered models that is expected to describe the phenomena 
of "high temperature superconductivity" in certain matterials |35],[36) • The motiviation for introducing the statistical 
gauge degrees of freedom there is quite different from the ones in this work. Here, the statistical gauge degrees 
of freedom are introduced a priori based upon locality and Lorentz invariance with the intention of describing the 
relativistic fermionic systems. Since the approach in this paper is applicable to all fermionic system, it is quite 
interesting to see whether the non-relativistic reduction of the problem can lead to some form of statistical gauge 
degrees of freedom for condensed matter systems at low energies, including the ones that Hubbard model describes. 
Nevertheless, many techniques in treating the statistical gauge degrees of freedom in condensed matter physics are 
expected to be applicable or a least adaptable here. 

One of the differences at formal level between the approach here and the ones used in condensed matter physics 
manifests in the different criterion for the selection of physical states within the full representing Hilbert space of 
the problem. The physical states in the statistical gauge theories developed so far in condensed matter physics is 
invariant under infinitesimal local gauge transformations, which is realized by the requirement that the operator form 
of the "Gauss law" annihilates all physical states in the superselection sector of the Hilbert space [37) . Such a strict 
enforcement of the statistical gauge invariance on the physical state vector of the system is neither necessary nor 
desirable for the statistical gauge invariant systems since it would exclude all finite density state from the physical 
sector of the system if one requires that the statistical "electric field" (denoted by ir u here) is finite. Such a situations is 
certainly unacceptable. Albeit this embarrassment can be circumvented in the 2+1 dimension | j37f , it is not expected to 
be easily implemented in higher dimensions. For the statistical gauge transformations, the statistical gauge invariance 
can be implemented by a imposing a less restrictive conditions on the physical superselection sector of the Hilbert 
space. Instead of requiring that the physical states are invariant under the gauge transformation, the gauge invariance 
on observables can be implemented by requiring that all states in a physical superselection sector of the Hilbert space 
cha nge a (coordinate dependent) common phase. The mathematical form for such a requirement is represented by 
Eq. 4.29. The physical superselection sectors arc then a functional determined by the common function c;. With such 
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a generalization of the "Gauss law" for the theory, both the requirements of the finiteness of the statistical "electric 
field" 7r u and of the fact that finite density states are actually physical states can be met consistently. 

The formalism is then applied to two half bosonized model Lagrangian densities. Four possible phases for the 
vacuum state of the interacting relativistic chiral symmetric systems are found. The first phase, called the O phase, 
correspond to the bare vacuum state of the system. Fermion-antifermion pairs condense in the second phase, named 
the a phase, of its vacuum. The a phase has the following properties: 1) the chiral SU(2)l x SU(2)r symmetry of 
the Lagrangian density of the system is spontaneously broken down to a SU(2)v symmetry 2) the baryon number 
density is zero 3) statistical blocking effects exists 4) the statistical gauge correlation is short ranged due to the 
presence of the statistical blocking effects. The third and fourth possible phases of the vacuum are called the ui phase 
and (3 phase respectively. Fermion pairs and antifermion pairs condense in the 10 and (3 phases. It is found that in 
these two phases of the vacuum: 1) the original chiral SU(2)l x SU(2)r symmetry of the Lagrangian density of the 
system remains unbroken in the to phase and is spontaneously broken down to a SU(2)y symmetry in the j3 phase 2) 
the baryon number density is different from zero or can be locally generated by separating fermion and antifermion 
rich region spontaneously 3) the U{1) symmetry corresponding to electromagnetism is spontaneously broken down 
to generate "massive photon" excitations (see ||) 4) no statistical blocking effects in these two phases 5) the spatial 
components of the statistical gauge excitation is long ranged; the quantum fluctuation in the time component of the 
primary statistical gauge field is Gaussian damped 6) off diagonal long range order exists in these two phases to give 
rise to macroscopic quantum behavior for the system, which is suppressed in the normal phase of the system. 

The implication of the finding presented in this paper on physical processes of strong interaction phenomena that 
are currently being or going to be observed or are in need to be explained theoretically remains to be investigated in 
the future. 
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FIG. 1. The p integration contour for the effective action. Here the filled circle represents the possible discrete spectra 
of S^lf] and the thick lines extending to positive and negative infinity represent the branch cuts of the logarithmic function. 
The ±iir above or under the thick lines denote the imaginary parts of the logarithmic function on the physical sheet of the p° 
plane. 

FIG. 2. Contour I is commonly used in literature where a non-covariant cutoff in 3-momentum space is applied. It is called 
the quasiparticle contour in this paper. Contour II is the one for the Euclidean effective action in the conventional approach. 

FIG. 3. The /i dependence of V e ff of the conventional approach in the a phase. It has two minima at non-zero /j,. Here A 
is the Euclidean momentum cutoff that defines the model. 

FIG. 4. The phase boundaries between the a, lu and the O phases. The chiral symmetry is unbroken in both the u> phase 
and the O phase. The a phase breaks the chiral symmetry spontaneously down to a flavor symmetry. 

FIG. 5. The phase boundaries between the a, (3 and O phases. The chiral symmetry is unbroken in the O phase. The a 
phase and j3 phase break the chiral symmetry spontaneously down to a flavor symmetry. 

FIG. 6. The p° integration contour for the effective action for transition amplitudes between states in which both the fermion 
and antifermion states with absolute value of their energy below e filled. Here the filled circle represents the possible discrete 
spectra of and the thick lines extending to positive and negative infinity represent the branch cuts of the logarithmic 

function. It also represents the p° integration contour of the full theory in its original Minkowski spacetime form. 

FIG. 7. The quasiparticle p° integration contour for the full theory is labeled by "I". The Euclidean p° integration contour 
for the full theory, which preserves the causal structure of the original one, is labeled by "II" . 

FIG. 8. The dependences of 14// (/it, e) in the a phase on /i and e along different directions in the fi-e plane. The direction 
in which V e // has the smallest value is in the /i = direction. Here v e ff = V e ///A 4 and x variable is either /i/A or e/A. 

FIG. 9. The dependences of V e //(/i, e) in the w phase on fi and e along different directions in the fi-e plane. The direction 
in which V e // has the smallest value is in the e = direction. Here « e // = T4///A 4 and x variable is either fi/A or e/A. 

FIG. 10. The dependences of V e //(/tt, e) in the /3 phase on fj, and e along different directions in the /i-e plane. The direction 
in which V e ff has the smallest value is in the e — direction. Here « e // = V e ///A 4 and x variable is either fi/A or e/A. 



1/3 

FIG. 11. The dependence of p vac on the spacetime independent part of the primary statistical gauge field /i. Here a = zuA. 
The unit for the dimensional quantities are GeV. Solid lines represent the case of free theory with mass and 0.5 respectively. 
Other lines represent the results for the massless NJL model with different strength of local fluctuations characterized by the 
a of the order parameter. 
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APPENDIX A: FREE FERMION SYSTEMS WITH n f FLAVORS AND n c COLORS 



1. With primary statistical gauge field only 



The Lagrangian density for a massive fermionic system with the primary statistical gauge held fx a (see Eq. 3.1) 
included has a form 

L = i*(^ + /i0 3 -m)tf, (Al) 

where m is the mass of the fermion. 

The generating functional W\fj, r], p] for such a system can be written as 

e w% v ,A = J D ^yI d H^+^) = cons t x e hSpi*n°iS™-*+fts<pi (A2) 
with rj and n Grassmann external helds and 



This equation allows us to write 



iS ( p ] 1 =ip + {iOz- m. (A3) 



iW[0, 0, fj] = ^SpLnrf°iS t p 1 + const, (A4) 



which for an uniform n a = (a, 0) , takes the following form 

d 4 p 



iW[0,0,fj]/Vi = -in trie / 77- 4t In (pi. - m 2 ) (pi - mr) + const, (A5) 
Jc ( 2n ) 

where V4 = L 3 T with L 3 — > 00 the spatial volume and T — > 00 the temporal extension of the system, p^_ = (p° + p, p), 
p^_ = (p° — fx, p) and C denotes the p° integration contour in the complex p° plane shown in Fig. ^. W[0, 0, fx] can be 
further specihed by requiring W[0, 0, 0] = 0. Such a VF[0, 0, is 

«m 0, an- / c 1^ - w ^ ) 7 m ' 1 - (a«) 

Since for free helds, the asymptotic grand canonic ensemble is the grand canonic ensemble and the quasiparticle 
approximation is actually an exact one (Appendix |b|) , the integration contour for p° integration above can be chosen 
as the one shown in Fig. ^ Q . Since the quantum fluctuations of the free fields are also exactly known and included 
already, they should not be sampled using the Euclidean approach. To avoid over counting of the quantum fluctuations 
of the free field, the p° integration should be done first and then the spatial component of p^. The adoption of the 
8-component spinor for the fermion held makes it equivalent whether the p° integration is carried out on the real axis 
or is on the imaginary axis || provided that it is done first. The result is 

iW[0,0,fA/Vi = ^ dpp 2 (u-E p ) = up-s = -n/V 4 , (A7) 

71" J0 



where e is the internal energy density and p is the fermion number density of the system and E p = \J p 2 + m 2 . It is 
evident that iW^O, 0, fi] correspond to the negative of the grand-potential at zero temperature in a many body 
system. 

Before ending of this subsection, it is worth mentioni ng t hat had we adopted the usual 4-component representation 



for a Dirac spinor, we would not have obtained Eq. A7 by using the path integration method with the primary 



statistical gauge held developed here. Part of the reasons is due to the non-symmetric way of introducing the primary 
statistical gauge field in the 4-component representation. This issue is discussed in in more details. 
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2. The full theory 



According to section [V, the generating functional W[r), to /j, e] for an uniform free fermion system in the full theory 
can be written as 



iW[rj, to fj,, e]/V4 = —inp 



f-q.p. ( 2?r ) 



d 4 p 

ln(p^ — m 2 )(p 2 _ — to 2 ) + const 



(A8) 



with "f.q.p." denoting the full quasiparticle p° integration contour in the complex p° plane shown in Fig. ^. The 
constant in the above equation is so chosen that W[0, 0, 0, 0] = 0, which means that 



iW[0,0, n,e]/V4 = —infUc 



d 4 p 



ln(pi -m 2 ){p 2 _ -m 2 ) -2 



d 4 p 



ln(p 2 — to 2 ) 



(A9) 



//.**. ( 2 ^) 4 

where "q.p." denoting the quasiparticle p° contour shown in Fig. 0. It is found, after some algebra, that 

iW[0, 0, /j, e]/Vi = (m+P(+) + M-P(-) - MPj - + £(-) - e) . ( Al 0) 



where /i± = /1 ± e and 



P(±) 
£(±) 



7T 2 



d|p|| P | 2 £; p 



(All) 
(A12) 



APPENDIX B: THE EXISTENCE OF THE DARK COMPONENTS FOR LOCAL OBSERVABLES 

1. Cluster decomposition and the origin of the dark component 

For simplicity, the two flavor half bosonized NJL model (see Appendix [e]) is used for our discussion. The model 
Lagrangian density with a primary statistical gauge field p a {x) is 

£ = 1^ (ifi + — o~ — i 7 5 f • 7?) * - -i- (a 2 + 7? 2 ) . (Bl) 
2 2Go 

After performing the path integration over the fermion fields ^ and 'J, the generating functional W[J] can be 
written as 

e w[J,ri = J D[a,if\e iS * ff[ '' M+i S dl 'f- J , (B2) 

where "J" and "/" represent, collectively, the external fields and auxiliary fields respectively. The effective action 
S e ff is given by 

1 1 f 

S e ff[a,TT,u] = -i-SpLnSp 1 [cr,Tf,u}S F [0, 0,0] + — - d 4 x (a 2 + n 2 ) , (B3) 
A ZUq J 

where Sp denotes the functional trace. The operator iSp 1 [u 1 tt, u] — ifi + (iO-$ — a — i7 5 r • 7? is the inversed propagator 
of the fermions in the background auxiliary fields o{x) and 7?(x). 

Let us define p[a, w , fj,; x] = Tr7°(cc |5f[ct, tt, u]\ x), where the trace "Tr" is over the internal degrees of freedom of 
the fermions. The fermion number density of the model in vacuum state is 

^mfJ](f,t = +™\m\f,t = -™) = \j^- ) =^ J D[<j,n}p[a,7T,n;x}e^ff^\ (B4) 

where Z = J D[a, jt}e lS ' lIf ^^'^ and / denotes the collection of {a, 7?} and W[/',/] is the weight functional of the 
asymptotic grand canonic ensemble discussed in the main text. The Minkowski spacetime is not suitable to study the 
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properties of the vacuum state using Eq. B4 since the initial and final auxiliary field configurations are not specified. 
The usual procedure to project out the contributions of the vacuum state is go to the Euclidean spacetime in which 
the vacuum state has lowest energy. 

In the mean field approximation, the vacuum phase is determined by minimizing the effective potential V e f / (a, /i) = 
—S e ff[a, 0, /Lt]/fi with Q — > oo the spacetime volume of the system, tr, /i Q spacetime independent and if assumed zero. 
Due to the chiral symmetry, the assumption 7? = does not result in a loss generality. The phase of the system 
is determined by the condition SV e ff(cr,iJ,)/6a — 0. The solution for a is non-zero after the coupling constant Go 
is greater than a critical value Gq c . A non-vanishing W generates an effective mass for the fermions, which act as 
quasiparticles. 



The mean field fermion number density for the vacuum is obtained from Eq. B4 by ignoring the functional 
integration over a and 7? and let a — a. The result is (see Appendix |Aj) 

Pmf = P{&, 0, n\ x) = -^-0(n - a) [n 2 - a 2 } , (B5) 

which is non-zero only when fi > a just like the fermion number density of free massive particles with mass m = a. 
Such a behavior of p is also predicted in the finite density field theory based on a global chemical potential /j, c h 
with N app discussed in the introduction behave in the same way as nj) MF . This will be discussed in the following. 
Therefore, the quasiparticle contributions saturate the fermion number density in the field theoretical approach to 
finite density problems based on a global chemical potential. 



The contributions of quantum fluctuations around the mean field a are formally included in Eq. |B4]. The results 
are commonly expressed as loop corrections to the fermion number density vertex, which is not attempted here. 

Instead of performing a loop expansion computation of the fermion number density, the effects of the quantum 
fluctuations can be evaluated non-perturbatively by "doing" the path integration. 

To proceed, the system under consideration is first putted in a Euclidean spacetime box of length L in each direction 
and with periodic boundary conditions at its boundary surfaces. The thermodynamic limit is defined as the limit 
of L — > 00. In the thermodynamic limit, the extremal configuration dominates the path integration among those 
configurations of a and 7? that give divergent action in the thermodynamic limit. The contributing finite action 
quantum fluctuation configurations, the number of which is proportional to the spacetime volume 11 = L 4 , are further 
classified into two categories: 1) correlated localized configurations, which are defined as the ones that approaches to 
the mean field configurations in the spacetime infinity and 2) correlated extended configurations, which are the ones 
that remain different from the mean field configurations at the spacetime infinity. 

For a given system, whether or not a configuration is a correlated extended configurations or is of localized ones^ is 
determined by dynamics. The on-shell amplitudes, are solutions of the "classical equation of motion" in the Euclidean 
spacetime 

-TFm = °' (B6) 
of{x) 

with S^ff the Euclidean effective action, / repres enti ng a or if fields. The superscript E shall be suppressed in the 



following. The set of the extended solutions to Eq. B6 are correlated ones. Albeit there are plenty of extended on-shell 
amplitudes in the Minkowski spacetime, there is no known one in the Euclidean one. I shall assume the absence of 
them. The degree of correlation of an arbitrary extended configuration at different spacetime points is determined by 
the degree of their deviation from the extended on-shell amplitudes for propagating excitations of the system. 

The correlation between the off-shell configurations at two different space time points decreases exponentially. For 
example, consider the field-field correlations or propagators 

(0\Tf( Xl )f(x 2 )\0) = / ^*{L e -*-(»i-« a >_L^ (B7) 

with T denoting time ordering. The off-shellness of these configurations is measured by their mass m. The corre- 
lation of these configurations at two different locations x\ and x 2 (in the Euclidean spacetime) decreases as fast as 
exp(— rm)/r with r — \x\ ^ x 2 \. So they can be decomposed into a superposition of localized ones with sizes of order 
I/to. For these set of configurations, one can divide the spacetime into cells with a dimension sufficiently larger than 



5 In the sense that it can be decomposed into localized ones. 
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their correlation length. Then the contribution of this set of field configurations to the partition functional Eq. B2 
can be cluster decomposed to 

Z[J]^l[z k [J], W[J] «5>*[J] (B8) 

k k 

with Wk[J] = lnzfe[J], Zk[J] and Wk[J] the corresponding partition functional of the kth cell. The full partition of the 
kth cell is 

z k [0] ~ / D k [f]e- S *^ 1+f \ (B9) 



where [/] is the effective action of the kth cell and the path integration of / is over those ones that equal to the 

mean field value / outside of the kth cell but with arbitrary amplitudes inside the finite volume. For a given theory, 
instead of arbitrary division of the spacetime, it is expected that there is an optimal one with minimum volume w 
for each cell and yet has an error below a predetermined one. We shall assume that such an optimal division of the 
spacetime into cells has already been found in the following discussion. 

Let us evaluate the contributions of the uncorrelated localized quantum fluctuations of the a and if fields to the 



vacuum fcrmion number density using Eqs. B4, B8 and B9 in the Euclidean spacetime. 



In the phase where the chiral SU(2)l x SU(2)r symmetry is spontaneously broken down, the "chiral angle" variable 
represented by 7? in the phase where (0 | -7? | 0) = becomes massless following the Goldstone theorem. The correlation 
length of the 7? field becomes divergent in the chiral symmetric limit. Therefore the field configurations of the 
Goldstone boson degrees of freedom contain the dominating on-shell components that can be included by doing a loop 
expansion as usual. Such a loop expansion contains no infrared divergences. Because the fermion number density 
p(a, 7? , p; x) under study is chiral symmetric, which means that it does not depend on a spacetime independent global 
"chiral angle"; it depends only on the derivatives of the "chiral angle" variables. The absence of a dependence of the 
fermion number density on a global "chiral angle" guarantees the absence of the infrared divergences in the quantum 
corrections from the Goldstone bosons. The configurations of the "chiral angle", being on the edge of their shell and 
extended configurations in nature, are uniform and infinitesimal in amplitudes since it contains no infrared divergences 
and has a number of distinct modes proportional to the volume £1 of the system. They can not modify the qualitative 
features of the quasiparticles. So, the 7? variable in the vacuum fermion number density can be eliminated. It is 
treated as zero in the following discussions. 

The quantum fluctuation in the "mass" term, namely the chiral radius or order parameter represented by a (when 
(0 \tt I 0) = 0) has different characteristics due to the fact that it contains no on-shell Euclidean configurations. These 
off-shell configurations have only short range correlations in spacetime. 

If the spacetime is divided into cells with their dimension optimally determined, then the path integration within 
each cell can be done independently. This gives us a cluster decomposed partition functional of the form given by Eq. 



B8 with the partition functional for each cell computed independent of each other. 

The cluster decomposition property of the pa rtit ion functional of the system reduces the full fermion number density 



of the vacuum given by Euclidean form of Eq. B4 to 

Pva~ 



~- J iW]/^W,0,M;z]e- S => + ^ ] , (B10) 
where the cell labeled by k is the one that contains the spacetime point x, J D^a 1 ] denotes integration over field 

(k) 

configurations that approach the mean field value outside the cell and S f f is the Euclidean effective action of the 
cell. " 



Eq. B10 is still too complicated to evaluate analytically. We make a further simplification by assuming that the 
functional integration of a' within a spacetime cell can be replaced by an ordinary integration over the spacetime 
averaged value of a' within that cell. It can be achieved by keeping the average of a' fixed while integrate over the 
rest degrees of freedom. The non-trivial part of it is in the assumption that after eliminating the rest of the degrees of 
freedom, the resulting effective potential V e ff remains, at least in form, the same as the original one. This procedure 



is in the same spirit as the renormalization group analysis. In this way, Eq. B10 reduces to 



1 r°° 

I d5ap[a + 5a,0,Mx]e-™ v °f^ +5 ^\ (Bll) 
z ["l J— 00 

where da =< a 1 > is the spacetime average of a' within the cell and the same reduction is also made to z[0]. The 
finiteness of w result in different qualitative behavior for p vac as a function of p. To explicitly see the difference, let 
us expand V e f / around a, keeping only the leading quadratic term 
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V ef f (a + 6a) = V ef f (a) + A(a)Sa 2 + ... (B12) 



and using Eq. B5 for the trace term. If wA is sufficiently large, the result can be written as 

f — ~o 



p vac = JLJ?A f _ d Sae~^ 2 [ M 2 - (a + 5af] 3/2 . (B13) 



It is clearly non-zero for any finite p below a since A is finite. The p dependence of p vac differs from the form given 



by Eq. B5 as long as Go is finite. Eq. B5 is recovered only in the limit of Go — > 0, which represents the free field case. 
The dependence of p\j a \ on p for a set of different values of a = vjA is plotted in Fig. O. Instead of a sharp rise in 
Pvac at p, = a, p is non-vanishing all the way to p, = 0. This component of the fcrmion number density can certainly 
not be attributed to the contributions of the quasiparticles. It is found therefore that there is a dark component for 
the fermion number density that can not be accounted for by the quasiparticle contributions. 

The NJL model has only one non-trivial vacuum phase since there is one independent absolute minimum in its 
effective potential. If the effective potential of the system contains a second local minimum with a higher energy 
density than the absolute minimum, then there are contributions from the local minimum to the fermion number 
density. The existence of such a contribution is another direct consequences of the cluster decomposability of the 
partition functional of the system. Models with a second minimum are studied in the literature. Those that have 
only one order parameter are represented by the Friedberg-Lee model |30[ ] in which the two minima of the effective 
potential correspond to confinement and deconfinement phase of the model. Those that has two order parameters are 
introduced in section |y|. 

In the presence of a higher virtual phase that is separated from the actual phase of the system by a potential barrier, 
the fermion number density of the system is saturated by both the quasiparticles of the first phase and those ones of 
the second virtual phase together with their corresponding dark component discussed above. It has a form 

= 1 + e - CT A » ( B14 ) 

where A = V^l — V^j > is the difference in energy density between the virtual phase and the actual phase, p vaCl 
and p V ac 2 are t ne vacuum fermion number density of the actual vacuum phase and that of the virtual vacuum phase 
respectively and vo is the optimal volume of the spacetime cell between which the order parameters of the system are 
discorrelated. The contributions of the virtual phase to the fermion number density or to any local physical observables 
are non-perturbative effects. Attempts had been made in Refs. |||l8],[34| to search for possible other virtual phases 
like the /3 or uj phase of the strong interaction vacuum. 



2. The reemergence of the quasiparticle picture 

The fermion number density discussed above are defined on a spacetime point. Such a precision is non- achievable 
in realistic observations. The physical observables can be represented by a "coarse-grained averaging" of the form 
~p vac = AN vac / AQ, with Ail the smallest volume in spacetime that the observation apparatus can resolve and AN vac 
the average number of fermions due to the coherent response of the system to an external classical field K within that 
volume. AN vac is not in general identical to J Af2 d A xp vac when Afl >> zu. It is obtained from the partition functional 
W[J, p] by adding the external field K to p. K has a constant non-zero value only within spacetime volume Af2, then 
AN vac — dW[0, p + K]/diK\K=a- In case when Ail << uj, namely the precision of the observation is much higher 
than the correlation length of the order parameter, the observed fermion number density ~p vac behaves in the same 
way as p vac - On the other hand, if Afl » w, then the smallest spacetime cell that contributes to AN vac is a region 
of volume Af2 rather than w. In this case, ~p vac is obtained from p vac by substituting AQ for vj. The effects of the 
dark component in the observed fermion number density are reduced as a result. When the resolution is sufficiently 
low, which means AHA >> I, the quasiparticle picture with only one vacuum , namely the actual one, reemerges in 
the response of the system to K. Therefore the dark component of the fermion number density is of transient nature. 
In most of the situations encountered in non-relativistic condensed matter system, the condition flA >> I is satisfied 
so that a global chemical potential approach is sufficient. 



3. The observation of the dark component at low energies 

On the other hand, the K field radiated by the vacuum fermion number density is of the form K(x) = 
J d 4 x'G(x,x')p vac (x') with G(x,x') the Green function for K (x). In a local QFT, it is p vac that is the source 
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for K(x) rather than p vac no matter how slow the resulting K(x) varies in spacetime. Thus the effects of the dark 
component are indirectly observable even in low energy processes. 

APPENDIX C: QCD LAGRANGIAN DENSITY WITH QUARK FIELD IN THE REAL 8 COMPONENT 

REPRESENTATION 

The QCD Lagrangian density for strong interaction can be written as 

C QCD = --TrG^G^ + ~W (ip - igA - m ) % (CI) 

where 

8 A a 

= d^A" - d v A tl - iglA^fA"], A» = ^B%— ^ C2 ^ 

and 



.,- 2 



A Q = \ [(1 + 3 )X a - (1 - 3 )A aT ] (C3) 

with mo current quark mass, g the QCD coupling constant, B a the bosonic gauge fields, "T" representing transpose 
and A a (a = 1, 2, . . . , 8) the Cell-Mann matrices. 

It is easy to verify that A a (a — 1, 2, . . . , 8) satisfy the same commutation relation with each other as the Gell-Mann 
matrices A a do, namely, 

[A Q ,A b ] =if abc A c , (C4) 

where f abc (with a,b,c — 1, 2, . . . , 8) are the set of group structure constants of SU(3). Therefore A a belong to the 
same adjoint representation of SU (3) as the one that A a belong. 

APPENDIX D: A FORMAL CONNECTION BETWEEN S eff [f] AND T[f] 

A formal approach to relate S e ff[f] and F[/] can be found using the method developed in Ref. jl4j by identifying 
S e ff[f] as I [</)]. Assuming that / are all real, which does not loss any generality since a complex field can be regarded 
as two real fields, the relation between T[f] and S e ff[f] can be established through a generalized vertex functional 
r[/, G], which is defined by 

f[f,G] = iS eff [f] - ^SpLnDG- 1 - ±Sp {V^G - l) + F 2 [f, G], (Dl) 

where T 2 [/, G] is the contributions of all two particle irreducible graphs with lines representing G in the shifted 
background fields / and vertices obtainable from S e ff[f] by expanding it around a set of shifted {/$} and keeping 
terms cubic in / or higher. Here D is the bare propagator for the boson fields and is symbolically given by 

v -l = $ 2 Seff{f] , , 

SfSf v ' 

The proper vertex function T[f] equals to T[f, G] with G satisfying the equation 

5T[f,G] 



5G 



= 0. (D3) 



In the Euclidean spacetime formulation, Eq. Dl becomes 



G] = S eff [f] - ^SpLnDG- 1 - ±Sp {V^G - l) + F 2 [f, G], (D4) 
where all the quantities above are evaluated in the Euclidean spacetime. 
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APPENDIX E: THE NAMBU JONA-LASINIO MODEL 



For the simplicity of the discussion, I consider a two flavor NJL model [ [l7| . In the conventional 4-dimensional 
representation for the fermion fields, it takes the following form 



C = ipiftip + Go (VV) 2 + (iprfrip) 



(El) 



with Go the coupling constant and f the three Pauli matrices in the isospin space. It has a SU(2)l x SU(2)r chiral 
symmetry. Due to the non-linear 4-fermion interaction term, it is not directly solvable. One of the best way to tackle 
the non- linear 4-fermion interaction model is to introduce auxiliary fields |8|. For this model two sets of auxiliary 
fields, namely a and tt are necessary. The model Lagrangian density after the introduction of the auxiliary fields a 
and 7? is of the following form 



1 — 



1 



£ = (ip - a - in ■ f 7 5 3 ) * - 777- K + tt 2 ) , 

which is written in terms of the 8-dimensional representation ^ for the Dirac spinors for fermions. 

The effective potential for this model can be computed using the Euclidean contour shown in Fig 
following form 



(E2) 



It has the 



47T 



1/1 6\ a 2 3 , / a 4 3, / a 3 



(E3) 



with a = G A 2 /Att. 

It has two phases. The system is in the first one, namely the O phase, when cxq < tt/12. It is in the second one 
called the a phase when the coupling constant ao > tt/12, where the original chiral SU(2)l x SU(2)r symmetry is 
spontaneously broken down to a SU(2)y flavor (isospin) symmetry and quark-antiquark pair condenses. Within the 
approximation adopted, the phase transition across the ao — tt/12 point is second order. 
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